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ABSTRACT 


A  brief  summary  of  the  important  results  obtained  under  Contract  AF-18 
(603) -100  is  given.  The  subject  matter  covered  by  the  investigations  included: 
(1)  high  explosive  generators  for  fast  particles,  (2)  mechanism  of  cratering 
in  ultra  high  velocity  impact,  (3)  observations  of  vaporization  accompanying 
ultra  high  velocity  impact,  (4)  mechanism  for  crater  expansion  in  shaped 
change  penetration,  (5)  ionization  and  electron  densities  in  detonating 
solid  explosives,  (6)  electrical  fields  and  electromagnetic  radiation  from 
chemical  detonations,  (7)  external  detonation  generated  plasmas,  (8)  the 
effect  of  pressure  on  the  degree  of  ionization  in  gaseous  detonations,  (9)  ioni¬ 
zation  and  electrical  conductivity  and  its  relationship  to  the  deflagration  to 
detonation  transition  in  solid  explosives,  and  (10)  transition  to  detonation  in 
liquid  explosives. 


INVESTIGATIONS  OF  HIGH  VELOCITY  IMPACT  AND  SOME  HIGH  EXPLOSIVES  PHENOMENA 


I .  Introduction 

The  following  constitutes  the  final  report  on  Contract  AF-18(603) -100  entitled 
"Investigation  of  High  Velocity  Impact  and  Some  High  Explosives  Phenomena".  Work 
commenced  on  the  contract  1  July  1955  and  continued  until  30  September  1961.  The 
Investigations  carried  out  during  this  period  may  be  considered  to  fall  under  the 
following  general  classifications:  (1)  high  explosive  generators  for  fast  particles, 
(2)  mechanism  of  cratering  In  ultra-high  velocity  Impact,  (3)  observations  of 
vaporization  accompanying"  ultra-high  velocity  Impact,  (4)  mechanism  for  crater 
expansion  In  shaped  charge  penetration,  (5)  Ionization  and  electron  densities  In 
detonating  solid  explosives,  (6)  electrical  fields  and  electromagnetic  radiation 
from  chemical  detonations,  (7)  external  detonation  generated  plasmas,  (8)  the 
effect  of  pressure  on  the  degree  of  Ionization  In  gaseous  detonations,  (9) 

Ionization  and  electrical  conductivity  and  Its  relationship  to  the  deflagration 
to  detonation  transition  In  solid  explosives,  and  (10)  transition  to  detonation 
in  liquid  explosives. 

The  following  constitutes  the  list  of  technical  notes  which  have  been 
published  as  a  result  of  research  on  these  subjects. 

Technical  reports  originating  from  work  under  Contract  AF- 18(603) -100: 

,1.  "Measurements  of  Ionization  and  Electron  Densities  in  the  Detonation  Waves 
of  Solid  Explosives",  Report  #1,  September  15,  1956. 

2.  "Behavior  of  Metals  under  Impact  Loading  and  the  Mechanism  of  Cratering", 
Report  #2,  October  15,  1956. 

3.  "Mechanism  of  Cratering  in  Ultra-High  Velocity  Impact",  10  July  1957, 
AFOSR-TN-57-486,  AD  136  479. 

4.  "The  Generation  of  High  Velocity  Projectiles  with  High  Explosives", 
AFOSR-TN-57-696,  AD  136  899,  1  October  1957. 

5.  "Propagation  Characteristics  of  Detonation  Generated  Plasmas",  AFOSR-TN- 
58-754,  AD  201  613,  20  June  1958. 

6.  "Mechanism  of  Cratering  in  Ultra-High  Velocity  Impact",  AFOSR-TN- 59-50 , 

AD  209  413,  30  January  1959.  I 

7.  "Electrical  Fields  and  Electromagnetic  Radiation  from  Chemical  Explosions", 
AFOSR-TN- 59 -551,  AD  216  691,  19  March  1959. 
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8.  Observations  of  Vaporization  Accon?>anying  Ultra-High  Velocity  Impact", 
AFOSR-ro-60-327,  4  January  1960. 

9.  "Framing  Camera  Observations  of  Ultra-High  Velocity  Penetration  in 
Transparent  Targets  and  a  Mechanism  of  Crater  Expansion",  AFOSR-TN-60-851, 

31  March  1960. 

10.  "Detonation  Generated  Plasmas",  AFOSR-705,  January  24,  1961. 

11.  The  Effect  of  Pressure  on  Ionization  in  Gaseous  Detonation",  AFOSR-704 

March  31,  1961.  ’ 

12.  "Ionization  Waves  from  Free  Surfaces  of  Detonating  Explosives",  AFOSR- 
1335,  September  7,  1961. 

Publications  based  upon  information  given  in  the  technical  reports  are  as 
follows: 

1.  "Generation  of  High  Velocity  Projectiles",  M.  A  Cook  and  R.  T.  Keyes, 

J.  Appl.  Phys.,  29,  1651  (1958). 

2.  "Mechanism  of  Cratering  in  Ultra-High  Velocity  Impact",  M.  A.  Cook, 

J.  Appl.  Phys.,  30,  No.  5,  725  (1959). 

3.  "Propagation  Characteristics  of  Detonation-Generated  Plasmas", 

M.  A.  Cook,  R.  T.  Keyes  and  L.  L.  Udy,  J.  Appl.  Phys.,  30,  No.  12,  1881  (1959). 

4.  "Detonation- Gene rated  Plasmas",  A.  Bauer,  M.  A.  Cook  and  R.  T.  Keyes, 

Proc.  Roy.  Soc.,  Vol.  A259,  508  (1961). 

5.  "Microsecond  Framing  Camera  Photographs  of  High  Velocity  Inpact", 

M.  A.  Cook  and  R.  T  Keyes,  Paper  presented  at  the  Third  Symposium  on  Hypervelocity 
Impact,  Chicago,  Illinois,  October  7-9,  1958. 

6.  "Framing  Camera  Observations  of  Ultra-High  Velocity  Penetration  in 
Transparent  Targets  and  a  Mechanism  of  Crater  Expansion",  R.  W.  Bartlett,  M.  A. 
Cook,  and  R.  T.  Keyes,  Paper  presented  at  the  Fourth  Symposium  on  Hypervelocity 
Impact,  Eglin  Air  Force  Base,  April  27-28,  1960. 

7.  "Ionization  Waves  from  Free  Surfaces  of  Detonating  Explosives",  A.  Bauer, 

M.  A.  Cook  and  L.  A.  Rogers,  Paper  presented  at  the  Western  States  Section  of  the 
Combustion  Institute,  Berkley,  California,  September  7-8,  1961. 

II .  Summary  of  Results 

A.  High  Explosive  Generators  for  Fast  Particles 

The  developments  under  this  phase  have  been  described  in  item  #4  of  the  list 
of  technical  notes.  The  method  which  evolved  employed  5  cm  diameter  Composition  B 
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charges  containing  wave  control  inserts,  Fig.  1  being  a  diagram  of  the  basic 
charge.  The  wave  control  insert  consisted  of  an  80°  apex  cone  of  height  2.5  cm  which 
was  composed  of  inert  material.  It  was  expected  that  the  most  suitable  material 


120’—  150’ 


Figure  1.  Basic  charge  used  in  the  wave  shaping  method 
of  accelerating  projectiles. 

from  which  to  fabricate  this  cone  would  be  glass.  However,  steel  proved  to  be 
slightly  superior,  probably  for  the  reason  that  the  cones  could  be  fabricated  with 
more  precision  when  made  of  steel.  The  precision  of  the  charge  components  was  found 
to  be  a  very  important  factor  affecting  the  performance.  The  function  of  the  wave 
control  insert  was  to  produce  an  inverted  detonation  wave  (concave-shaped  looking 
from  the  front)  rather  than  the  usual  spherical  shaped  wave  (convex  toward  the 
front).  The  above  type  wave  control  insert  in  contrast  to  more  bulky  "interrupter 
type"  wave  shapers  which  shape  the  wave  by  stopping  detonations  completely  beyond 
the  insert,  after  which  the  wave  is  reformed  by  detonation  of  the  explosive  around 
the  sides  of  the  insert,  permits  initiation  of  detonation  by  the  shock  propagating 
through  the  insert.  With  such  a  device  the  variable  delay  desired  across  the  wave 
front  was  produced  by  a  combination  of  two  factors.  Firstly,  because  the  detonation 
velocity  of  Conposition  B  exceeded  the  velocity  of  the  compressional  wave  generated 
in  wave  shaper,  the  wave  was  delayed  progressively  more  toward  tha  axis  as  the 
path  length  through  the  shaper  increased.  The  second  contribution  resulted  from 
the  variable  delaying  of  initiation  of  the  explosive  in  the  wafer  beyond  the  wave 
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control  insert,  the  delay  of  initiation  increasing  as  the  path  length  through  the 
wave  shaper  increased  which  resulted  because  the  attenuation  of  the  wave  was  greater. 
When  wave  shaping  is  accomplished  in  this  manner  the  impulse  of  the  detonation  wave 
does  not  appear  to  be  degraded  to  any  great  extent  as  in  the  case  for  "interrupter 
type"  wave  shapers  where  the  detonation  wave  has  to  reform  completely. 

The  benefit  of  such  an  "inverted"  wave  for  accelerating  projectiles  was 
believed  to  be  twofold.  Firstly,  an  inverted  detonation  wave  is  unstable.  The 
wave  rapidly  reverts  to  the  normal  spherical,  convex  shape  by  virtue  of  an  increased 
velocity  along  the  central  part  of  the  wave,  and  consequently  the  pellet  is  acceler¬ 
ated  by  a  wave  traveling  at  a  velocity  higher  than  the  normal  detonation  velocity. 
Secondly,  an  Inverted  wave  was  found  to  lead  to  reduced  fragmentation  of  the  pellet. 

Tests  were  carried  out  firstly  to  determine  the  most  suitable  shaped  pellet, 
l.e.,  the  shape  primarily  which  possessed  the  smallest  tendency  to  fragment.  The 
best  combination  was  found  to  be  a  120°  -  150°  pellet  (interior  angle  rear  surface 
120°  -  front  surface  150°  concave  inward)  or  Pugh  type  fragment.  Nextly,  using 
0.95  g,  1.9  cm  diameter  aluminum  pellets  of  this  type  and  80°  2.5  cm  height  conical 
shaped  wave  control  inserts  the  wafer  length  x  (and  consequently  the  shape  of  the 
wave  impacting  the  pellet)  was  varied  to  determine  the  influence  upon  projectile 
velocity  and  projectile  fragmentation.  Fortunately,  the  tendency  for  projectiles  to 
fragment  appeared  to  be  minimum  for  the  wafer  length  that  led  to  the  highest 
velocities  which  was  about  4.5  cm  using  either  glass  or  steel  wave  control  inerts. 
Under  these  optimum  conditions  velocities  up  to  7600  m/sec  were  realized  with  the 
0.95  g  aluminum  pe llets . 

B-  Mgchanism  of  Cratering  in  Ultra-High  Velocity  Impart 

This  subject  was  treated  in  items  #3  and  H  of  the  list  of  technical  notes. 
Briefly,  cratering  was  classified  into  three  regions  according  to  the  striking 
velocity,  as  follows: 

a.  Impact  at  a  striking  velocity  V  in  the  range  V  >  V>  C  where  C  is  sonic 
velocity  in  the  target  and  is  the  threshold  velocity” required  to  cause  plastic 

deformation,  i.e.,  the  velocity  at  which  the  Bernoulli  pressure,  1/2  p  (V-U)^  is 
just  equal  to  the  yield  strength  a,  being  the  density  of  the  projecUle  and 
U  the  penetration  velocity.  This  condition  is  possible  only  when  V  >  C  which 
requires  generally  a  relatively  low  acoustic  velocity  such  as  a  loose  sand  target. 


b.  Impact  at  a  striking  velocity,  V,  such  that  <V  <V^j.  where  is  the 
minimum  velocity  required  for  an  impact  explosion  of  the  target,  i.e.,  vaporization 
of  the  target.  The  type  (b)  penetration  (with  <  C)  is  probably  best  known 
from  studies  of  penetration  of  targets  by  shaped  charge  jets.  This  type  cratering 
has  been  treated  by  applying  the  laws  of  hydrodynamics  and  considering  the  collision 

process  under  the  approximation  of  non-viscous  fluid  flow. 

c.  Cratering  with  V  >V^j..  Impact  explosions  were  expected  to  occur  only 

for  striking  velocities  somewhat  above  and  C. 

The  type  (a)  penetration  was  not  treated  further  in  our  studies.  However,  the 
idealized  theory  of  type  (b)  penetration  was  reviewed,  and  then  the  principle  of 
conservation  of  momentum  was  applied  in  a  new  way,  yielding  in  the  case  of  pene¬ 
tration  by  shaped  charge  jets  an  expression  for  the  final  crater  diameter.  The 
reasoning  was  as  follows;  During  an  infinitesmal  increment  of  penetration  &P 
which  occurs  during  the  time  6t  a  total  impulse 

61  =  l/2XPj(V^-U^)\6t  (1) 

is  imparted  to  the  target  where p ^  is  the  density  of  the  jet,  is  the  jet 
velocity,  U  is  the  penetration  velocity,  is  the  cross-sectional  area  of  the 
jet,  and  X  is  a  factor  to  account  for  gaps  between  particles  in  a  discontinuous 
jet.  (For  a  continuous  jet  x=  1.)  The  relation  1/2  Xp  j(V^-U^)^  is  simply  the 
Bernoulli  pressure  generated  by  a  jet  of  velocity  in  fluid  flow  penetrating  a 

target  at  velocity  U  .  Except  at  the  surface  of  the  target  where  relief  of  pressure 
allows  some  backward  flow,  this  increment  of  impulse  produces  effectively  only 
lateral  residual  flow.  Since  the  area  resisting  the  pressure  of  the  lateral 
residual  flow  at  the  collision  interface  according  to  this  concept  is  A^,  it 
was  natural  to  conclude  that  the  area  resisting  the  lateral  residual  flow  at  a 
particular  site  may  be  taken  to  be  A,  the  expanding  cross-sectional  area  of  the 
crater  at  the  particular  location  where  the  increment  of  impulse  was  first 
received.  Thus,  ^ 

61  =  1/2XP  j  ^o  °  ^t^ 

where  the  zero  subscripts  indicate  initial  conditions.  The  expression 
(1/2  p  +  (j)  is  simply  the  pressure  of  lateral  residual  flow,  p  ^  being  the 
target  density,  a  the  dynamic  yield  strength  of  the  target  and  U  the  lateral  flow 
velocity  at  the  crater  location  whose  instantaneous  cross-sectional  area  is  A. 
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The  crater  will  continue  to  expand  until  U  =  0,  i.e.,  until  the  pressure  of  lateral 


residual  flow  equals  the  static  yield  strength  of  the  target, 


Setting  U  *  0, 


o',  eliminating  U^,  and  expressing  cross-sectional  areas  of  the  crater  and  jet 


in  terms  of  their  diameters  led  to  the  expression 

d  p  (Xp,)l/2 

j  I  _  o  t  .1 


(3) 


for  the  final  crater  diameter  in  terms  of  the  jet  diameter,  d  . 

An  expression  for  crater  volume  in  single  particle  type  (b)  impact  was  developed 
where  the  energy  radiated  from  the  collision  zone  by  shock  waves  was  neglected. 

The  development  was  based  upon  the  idealized  situation  undeV  which  a  target  was 
considered  to  be  penetrated  by  a  cylindrical  projectile  of  length  and  diameter 
dp.  Under  these  conditions  a  cylindrical  crater  of  constant  radius  would  be 


formed,  the  radius  of  which  could  be  calculated  by  means  of  equation  3 
the  depth  of  penetration  being  given  by 

''  - 

which  is  a  well-known  equation  from  shaped  charge  theory, 
then  derived  to  be 

1  /•>  ^  n  2 

m  V 

1/2.2  20 


With 


o  1/2  ^  1/2 

% _ _ 


<P  +  Pt 


(4) 

The  crater  volume  was 


(5) 


■) 


which  for  p  =  reduced  to 
P  t 


m  V 
8o 


(6) 


where  m^  is  the  projectile  mass,  V  its  velocity,  and  o  the  yield  strength  of  the 
target.  For  a  spherical  particle  the  impulse  flow  of  the  target  would  be  radial 
and  the  resulting  crater  nearly  hemispherical.  However,  it  was  concluded  that 
the  geometry  of  the  dynamic  pressure  to  static  pressure  conversion  should  be  such 
that  equations  (5)  and  (6)  should  apply  for  spherical  particles  also. 

Type  (c)  cratering  was  then  discussed  and  predictions  were  made  concerning 
the  threshold  velocity  required  for  impact  explosions.  It  was  considered  that 
the  impact  velocity  would  be  such  a  velocity  that  the  heat  developed  by 
inelastic  collision  in  the  target  equaled  the  heat  of  vaporization,  and  that  in 
general  the  penetration  velocity  U  would  be  greater  than  the  shock  velocity 
corresponding  to  the  Bernoulli  pressure  at  the  collision  interface,  i.e.,  the  shock 
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wave  would  not  move  away  from  the  collision  interface,  and  thus  there  would  be  no 
means  of  removing  from  the  target  the  heat  of  inelastic  impact. 

Through  application  of  a  theory  for  the  compressibility  of  solids  the 
specific  compressional  work  of  inelastic  impact  was  calculated  in  terms  of  the 
impact  velocity  and  equated  to  the  heat  of  vaporization  to  yield  an  expression 
for  V  ,  the  threshold  velocity  for  vaporization,  the  result  being 

rou  -1/^  1/2  ^  .  1/2 

rsh  p-  p. 


which  for  equal  target  and  projectile  materials  reduces  to 

V  .  2(8h  /P  P)''‘  (8) 

ct  CO 

where  h  is  the  spacific  heat  of  vaporization  and  is  the  compressibility  at 
atmospheric  pressure  Some  typical  values  of  for  various  projectile-target 

combinations  calculated  by  means  of  these  equations  are  listed  in  Table  1. 

Loss  of  energy  from  the  material  of  the  target,  which  is  at  any  instant  in 
collision  with  the  projectile  by  means  of  the  shock  wave,  was  considered  to  occur 
only  when  U  <  where  U  is  the  penetration  velocity  and  is  the  velocity  of 

the  shock  wave  in  the  target  T»-,erefore,  the  radiated  shock  waves  should  tend 
to  prevent  an  impact  explosion  as  long  as  the  condition  U  <  maintains.  This 

situation  was  not  considered  to  apply  to  projectiles,  however,  because  shock 
waves  propagating  back  into  the  projectile  from  the  collision  zone  would  not  remove 
the  inelastic  impact  energy  from  the  projectile.  On  the  basis  of  the  assumption 
that  the  energy  accumulated  in  the  projectile  part  of  the  collision  zone  plus 
that  which  flowed  away  from  the  collision  zone  into  the  projectile  was  approximately 
the  same  as  the  energy  accumulating  in  the  target  material  plus  that  which  flowed 
into' the  target  in  radiated  shock  waves, calculations  showed  that  projectile  impact 
explosions  could  possibly  occur  at  velocities  only  about  45%  of  those  required  to 
produce  target  impact  explosions 

An  expression  was  next  derived  for  the  maximum  amount  of  material  which 
would  be  vaporized  from  a  target.  In  the  case  of  equal  target  and  projectile 
materials  the  relationship  in  terms  of  the  ratio  of  the  mass  of  target  material 
vaporized  to  the  projectile  mass  was  2 

'"t^’"p  -  It  ■  ^ 
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Table  I:  Critical  impact  and  penetration  velocities  for  target  impact  explosions 


Projectile 


Steel 

Aluminum 

Copper 

Magnesium 

Zinc 

Lead 

Steel 

II 

II 

II 

II 

Aluminum 

II 

It 

It 

II 

Copper 

It 

M 

II 

It 

Magnesium 

II 

It 

It 

It 

Zinc 

ti 

It 

It 

tt 

Lead 

It 

It 

tl 

tl 


Target 

h 

c 

P  o 

^ct 

ct 

^ct 

Kcal/g  atm"  ^*10 

megabars 

km/sec 

km/ sec 

Steel 

1.7 

5.94 

1.37 

11.8 

5.9 

Aluminum 

2.04 

13.85 

0.58 

13.0 

6.5 

Copper 

Magnesium 

1.28 

1.50 

7.27 

29.5 

1.14 

0.27 

10  0 
11.2 

5  1 
5.6 

'3  C 

Zinc 

0.41 

13.0 

0.43 

7 .0 

J  .  J 

Lead 

0.23 

24.2 

0.30 

4 . 6 

2.3 

Aluminum 

0.58 

10.3 

6.5 

Copper 

Magnesium 

Zinc 

1.14 

0.27 

0.43 

10.4 

8.2 

6.8 

5  1 
5.6 
3.5 

Lead 

0,30 

5.1 

2 . 3 

Steel 

1.37 

16.0 

5  .9 

Copper 

Magnesium 

Zinc 

1,14 

0.27 

0,43 

14.4 

10.1 

9.2 

5  1 
5.6 
3.5 

Lead 

0.30 

7.0 

2.3 

Steel 

■  1.37 

11.4 

5.9 

Aluminum 

0.58 

18.0 

6.5 

Magnesium 

Zinc 

0.27 

0.43 

8  •  1 
6.6 

5 . 6 
3.5 

Lead 

0  30 

4  9 

2.3 

Steel 

1.37 

18.5 

5.9 

Aluminum 

0  58 

14.6 

6.5 

Copper 

Zinc 

1.14 

0  43 

16.6 

10.6 

5.1 

3.5 

Lead 

0  30 

8.2 

2.3 

Steel 

1.37 

12.1 

5.9 

Aluminum 

0.58 

10.5 

6.5 

Copper 

Magnesium 

Lead 

1.14 

0.27 

0.30 

10.8 

8.4 

5.2 

5.1 

5.6 

2.3 

Steel 

1.37 

10  8 

5.9 

Aluminum 

0.58 

9.7 

6.  5 

Copper 

Magnesium 

Zinc 

1.14 

0.27 

0.43 

9  6 
7.8 
6  3 

5.1 

5.6 

3.5 

of  sound  - 

5.3  (steel) 

;  5.1  Al); 

3.6  (Cu); 

3.7  (Zn); 

1.23  (Pb) 

C  =  velocity  of  sound 
4.6  km/ sec  (Mg) . 
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The  second  paper  concerned  with  the  mechanism  of  cratering  in  ultra-high 

velocity  Impact,  l.e..  Item  #6  on  the  list  of  technical  notes,  was  concerned 

mainly  with  theoretical  evaluations  of  dissipative  factors  occurring  In  the 

cratering  process  and  the  effect  upon  the  Idealized  theory  discussed  In  the 

preceding  paper.  In  the  non-ldeallzed  case  under  the  conditions  of  plastic  flow 

continuity  of  pressure  at  the  Impact  junction  between  the  projectile  and  the 

target  may  be  expressed  by  the  equation 

Pp  (V-U)^/2  -  Pj.U^/2  +  ♦(U,a)  (10) 

where  ^(U,o)  Is  a  function  describing  dissipative  pressures,  e.g.,  losses 

associated  with  yield  strength,  compression,  heating  of  the  target,  radiating 

shock  waves,  and  possibly  Impact  explosions.  The  term  ^  was  therefore  the  dynamic 

yield  strength  of  the  target.  The  penetration  equations  for  shaped  charge  jets 

were  then  derived  for  the  non-ldeallzed  case,  namely 

2  1/2 

U  -  V  {  1  -  (A  +  2Vpj  (1-A)]^'  /(1t^))  (11) 


L  {  l-(  A+  2«/  y  (l-A)l^^^ 


(A+  2^/  (1-A)] -  A 


U  being  the  penetration  velocity,  V  the  jet  velocity,  P  the  depth  of  penetration, 
and  L  the  length  of  the  jet  (considered  to  be  of  uniform  velocity).  And  In  the 
non-ideal  case  the  equation  for  crater  volume,  equation  6,  became 

T-  T/4?  (13) 

where  T  is  the  kinetic  energy  of  the  projectile  just  before  it  strikes  the  target. 

In  most  of  the  treatments  of  penetration  by  shaped  charges,  except  that  of 
(2) 

Eichelberger,  who  took  into  account  the  static  yield  strength  of  the  target, 
has  been  considered  negligible.  The  study  made  at  this  laboratory  pointed  out 
that  the  Idealized  theory  has  given  good  results  when  one  Is  well  within  the 
realm  of  type  (b)  penetration.  This  theory,  however,  was  expected  to  break  do\m 
in  the  realm  of  impact  explosions  on  the  high  velocity  side  and  on  the  other  hand 
also  when  the  velocity  was  reduced  to  approach  the  threshold  of  plastic  deformation. 

Theoretical  developments  of  the  ♦  and  ♦  functions  were  next  carried  out.  The 
dynamic  yield  strength,  0  ,  was  considered  to  consist  of  four  terms  as  follows: 

4>*h+g+S  +  a  (14) 


i. 
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where  h  is  the  specific  adiabatic  work  of  compression  on  the  target  material, 
g  is  the  permanent  specific  shock  heating  of  the  target  material  resulting  from 
the  difference  in  entropy  between  a  shock  compression  and  an  adiabatic  compression, 

S  is  the  energy  radiated  from  the  collision  zone  by  shock  waves,  and  o  is  the  static 
yield  strength  of  the  target.  Through  application  of  an  equation  of  state  for 

(3) 

metals  the  h  and  g  functions  were  evaluated  to  be  as  follows: 

h  =  ap^  ^  (x^/a  -  2.3  log  x  -  a  (15a) 

g  -  k  (x  -  l)^/x^  (15b) 

where  p  is  the  compressibility  at  atmospheric  pressure,  x  =p  /p^  ,  a  is  a  slowly 
o  ^  o 

varying  parameter  which  for  most  metals  has  values  which  range  between  2  and  2.15, 

^  A  A  A 

k  is  a  constant  computed  to  be  1.88  x  10  ,  1.53  x  10  ,  1.84  x  10  and  1.26  x  10 

°K  for  iron,  lead,  copper,  and  aluminum,  respectively,  and  is  the  average  heat 

capacity  per  unit  volume  of  the  target  material.  The  density  ratio  before  and  after 

compression  was  related  to  the  pressure  through  the  relations 

P  /Po  =  (Pt^/t^t^^  (16a) 

P  =  (p^/p‘^/ap^  (16b) 

It  should  be  stressed  here  that  the  effective  pressure  used  in  computing  the 

2 

density  ratio  was  the  Bernoulli  pressure  of  flow,  P^U  /2. 

No  way  of  evaluating  S  was  worked  out.  However,  comparisons  between 
measured  and  theoretical  (based  upon  the  non- Idealized  theory)  U/V  vs  V  curves  showed 
that  S  was  only  important  in  very  low  velocity  impact,  i.e.,  at  about  the  threshold 
for  plastic  deformation  of  the  target.  One  would  expect  this  result  also  from 
framing  camera  pictures  taken  of  shaped  charge  jets  penetrating  transparent  targets 
which  showed  that  in  supersonic  penetration  the  shock  wave  did  not  propagate  into 
the  target  ahead  of  the  collision  zone  and  thus  did  not  radiate  energy  from  the 
collision  zone.  This  occurred  only  for  subsonic  penetration.  The  static  yield 
strength,  o  and  h,  the  adiabatic  work  of  compression,  were  found  to  be  the  largest 
contributors  to  <>.  In  the  case  of  steel,  for  example,  g  is  generally  about  0.1 
as  great  as  h 

More  exact  calculations  for  the  threshold  velocity,  to  produce  target 

explosions  were  performed,  being  taken  to  be  the  velocity  where  the  internal 

heating  of  the  target  by  adiabatic  work  of  compression  plus  shock  heating  equaled 
the  heat  of  vaporization  of  the  metal,  i.e.. 
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hc  =  (ap^)‘^t  x^/a  -  2.3  log  x  -  a'^]  +  k(x-l)^/x^  (17) 

Some  representative  results  are  given  in  Table  II  which  may  be  compared  with  those 
of  Table  1. 


Table  11- 

Threshold  conditions  for 
theory 

in^jact  explosions 

according 

to  non-idealized 

Metal 

h 

^ ,  2 
dynes/cm 

X  10*“ 

P  /p 
o  c 

^ct  2 
dynes/cm 

X  10-11 

^'ct 

km/ sec 

(km/sec) 

(steel 

projectile) 

Iron 

5  2 

0  ,56 

23  8 

6  9 

14.7 

Copper 

4  7 

0  55 

21  8 

6  2 

13.7 

Lead 

1.1 

0.57 

5  5 

2.8 

6.6 

Zinc 

1.3 

0  62 

4.7 

3.1 

6.5 

Aluminum 

2.3 

0  o6 

10  9 

8  0 

13.3 

Magnesium 

1  1 

0  58 

5.6 

7.2 

11.0 

Note  that  as  far  as  p  ^  and  U  ^  are  concerned,  the  results  in  Tables  I  and  II 

ct  ct  ’ 

were  not  exceedingly  different  However,  was  significantly  different. 

Crater  volume  in  ultra-high  velocity  impact  (type  c  cratering)  was  then 
treated  theoretically  Under  this  type  cratering  the  crater  was  considered  to 
be  divided  into  two  zones,  a  zone  (a)  which  comprised  the  part  of  the  target  where 
vaporization  or  impact  explosion  occurred,  and  zone  (b)  which  comprised  the  region 
of  plastic  deformation  The  crater  made  in  zone  (b)  was  considered  that  which 
would  be  created  by  an  explosive  detonated  in  the  volume  whose  available  energy 
was  T  ,  where  V  and  T  were  related  through  the  expression 

\  -  \\  -  V,0,u\,/2  (18) 

V  h  being  the  energy  lost  in  vaporization  of  zone  (a)  of  the  target,  and 

a  c  2 

V  pU  /2  being  the  energy  available  for  secondary  penetration  in  zone  (b) .  (Note 

3i  L 

that  equation  (13)  neglects  heat  loss  in  the  projectile.)  The  volume  in  zone  (b) 

was  then  taken  to  be  that  given  according  to  equation  (13),  i  e., 

=  T  /Ci 
b  a 

Solving  for  from  equation  (18)  yielded  the  equation 

V  =  mV^/2  (h  +  p  U  ^/2)  =  mV^/2n 
a  c  t  c  t  t^c 


(19) 

(20) 
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The  total  crater  volume  was  then  determined  to  be  ^ 

-  Va  +  Vt  .  Va  (1  +  Pt  ^et  “  ''a^t^ct 

The  volume  of  zone  (a)  was  shown  to  be  negligible  compared  to  that  of  zone 


(21) 

(b) 


2  — 

even  when  V  »  since  »  Sa¬ 

lable  III  lists  some  representative  ♦  /p  vs  p  data  for  steel  penetrating 
mild  steel  targets  which  was  calculated  using  equation  (14)  (neglecting  S)  and 


equations  (15a)  and  (15b). 


Table 

n 

111: 

Computed  ♦/p 

4) 

vs  p  data  for 

4>/p 

steel  penetrating 

P 

mild  steel  targets.* 

4>  O  /P 

r 

15.2 

5.2 

0.34 

0.6 

0.13 

0.22 

10.0 

3.0 

0.30 

0.4 

0.092 

0.23 

8.0 

2.3 

0.29 

0.2 

0.05 

0.25 

6.0 

1.6 

0.27 

0.1 

0.03 

0.30 

4.0 

1.0 

0.25 

0.08 

0.028 

0.35 

2.0 

0.48 

0.24 

0.06 

0.026 

0.43 

1.0 

0.23 

0.23 

0.04 

0.023 

0.60 

0.8 

0.18 

0.23 

,  2 

0.02 

0.02 

1.0 

2 

*  and  p  expressed  in  dynes/cm  *10 


One  will  note  that  Wp  varied  between  the  relatively  narrow  limits  of  0.34  to 
0  22  from  p  .to  0.04  p  .  Thus,  it  was  concluded  that  for  type  (b)  cratering  30  to 
40  percent  of  the  primary  impact  energy.  p//2,  would  be  lost  as  heat  in  residual 
flow  and  the  rest  would  be  lost  in  overcoming  the  yield  strength,  o  .  On  this  basis 
iT/owas  concluded  to  be  roughly  1.4-1. 7.  This  approximated  the  experimental  value 
of  */o  -  1.75  which  had  been  determined  for  steel  jets  penetrating  steel  targets. 
Thus,  it  was  concluded  that  the  approximation  ♦  -  1.75 o could  be  used  in  the 
equation  -r  -  T/4'5  for  the  crater  volume  in  type  (b)  cratering. 


C  observationn  of  Vaporization  Accompanying  Ultra-High  Velocity  Inyac_t 
The  results  of  this  study  were  described  in  item  #8  of  the  list  of  technical 
notes.  Briefly,  experiments  were  carried  out  at  velocities  up  to  8.5  km/sec  using 
shaped  charge  jets  and  explosively  propelled  projectiles  with  a  variety  of  target 
and  projectile  materials.  Initially,  the  impacts  were  recorded  by  microsecond 
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color  framing  camera  photographs.  These  photographs  revealed  that  at  sufficiently 
high  impact  velocities,  depending  upon  the  target  and  projectile  materials  Involved 
vapors  uere  generated  at  Impact.  These  vapors  were  highly  luminous,  their  luminosity 
being  found  to  result  primarily  from  the  decay  to  lower  energy  states  of  the  exited 
and  Ionised  gaseous  metallic  species  making  up  the  vapor,  rather  than  chemical 
reaction  with  the  atmosphere  surrounding  the  Impact  tone.  This  conclusion,  which 
proved  that  the  clouds  were  a  vapor  rather  than  fine  solid  particles,  was  based 
upon  the  fact  that  (1,  Impact  generated  vapors  were  found  through  probe  conducUon 
measurements  to  he  highly  conducting  and  therefore  Ionised  to  an  appreciable  extent. 

onseguently,  the  term  Impact-generated  plasmas  was  applied  to  these  luminous.  Impact¬ 
generated  Clouds.  (2,  ^e  l„lnoslty  of  these  l^.ct-generated  clouds  was  l^ger 
lasting  and  apparently  Just  as  Intense  when  the  Impact  took  place  In  as  near  perfect 

“  -“»^h»ere  such  as  oxygen. 

wh  ch  should  react  chemically  with  the  material  of  the  cloud.  In  general,  were  no 

rlghter  than  Impacts  In  relatively  Inert  atmospheres.  And  (4)  Impacts  Involving 

materials  possessing  low  cohesive  energies  such  as  lead  (which  one  would  expect  to 

vaporise  most  easily,  produced  more  Intense  plastm.  than  materials  which  should 

e  more  difficult  to  vaporUe  hut  which  would  tend  to  react  more  energetically  with 
the  atmosphere. 

Techniques  also  were  applied  which  collected  small  quantities  of  condensed 
vapors,  the  arrangement  being  such  that  the  likelihood  of  collecting  solid  particles 
was  mlntmlaed.  Such  samples  were  analysed  quantitatively  to  determine  the  ratios 
p  ojectlle  or  Jet  to  target  metals  In  the  vapor  condensate  as  a  function 
o  ^act  velocity.  In  most  cases  the  material  collected  consisted  of  metal  oxides 
hut  the  X-ray  diffraction  patterns  did  not  permit  conclusive  Identification  of 
specific  ccxnpounds.  These  studies  showed  that  for  all  the  systems  Investigated  the 
mete  w  t  the  lower  heat  of  vaporisation  was  more  abundantly  vaporised,  a  tact 
which  agreed  with  the  results  of  the  photographic  studies.  Studies  with  copper 
and  steel  used  alternately  for  the  Jet  and  the  target  Indicated  that  as  tar  as 
vaporisation  was  concerned.  It  did  not  matter  whether  a  particular  «tal  made 
up  the  target  or  the  Jet.  In  all  cases  more  copper  vaporised  than  Iron.  This  fact 
was  not  in  agreement  with  the  prediction  that  projectiles  would  vaporise  more  easily 
than  targets.  Another  dlagreement  with  the  predictions  of  the  theory  descrlhed 


i 
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earlier  pertained  to  the  threshold  for  vaporization.  In  practice  vaporization  was 
found  to  occur  at  lower  velocities  than  calculated,  although  this  was  not  surprising. 

In  instances  where  substitutional  solid  solution  alloys  were  possible  such 
alloys  could  form  in  the  impact  region.  For  example,  a  copper  target  was  Impacted 
by  a  steel  jet  after  piercing  a  thin  3%  zinc-aluminum  alloy  sheet  Inserted  in  front 
of  the  target.  The  surface  of  the  copper  target  was  shown  by  X-ray  fluorescence  to 
contain  zinc.  Furthermore,  the  most  intense  diffraction  peak  (111)  for  copper  was 
shifted  through  a  double  Bragg  angle  of  0.47°.  The  peak  was  .less  sharp  than  for 
pure  copper  and  was  also  asymmetric  indicating  that  the  lattice  expansion  was  not 
uniform  throughout.  Nevertheless,  if  as  a  first  approximation  the  location- of  the 
peak  maximum  was  assumed  to  represent  the  average  lattice,  an  average  lattice  expansion 
of  0.039  X  was  obtained  which  is  equivalent,  to  a  substitution  of  67.  zinc  in  the 
copper  lattice  near  the  surface.  Thus,  simultaneous  condensation  of  zinc  and  copper 
vapor  on  the  surface  of  the  copper  target  resulted  in  the  formation  of  brass. 

It  was  concluded  that  a  more  reliable  method  of  sampling  the  species  from  the 
impact  generated  plasmas  was  required  before  the  mechanism  governing  the  vaporization 
under  hypervelocity  impact  could  be  determined  experimentally.  However,  the  study 
did  serve  to  show  that  significant  amounts  of  vaporization  did  occur  under  proper 
conditions. 


D-  Mechanism  for  Crater  Expansion  in  Shaped  Charge  Penetration 
The  results  of  this  study  were  given  in  item  #9  of  the  list  of  technical  notes. 
This  study  consisted  firstly  of  extending  the  model  developed  in  item  #3  of  the 
list  of  technical  notes  which  predicted  final  crater  diameters  resulting  from  the 
penetration  of  solids  by  shaped  charge  jets,  the  extension  being  the  development  of 
an  equation  which  predicted  the  radius  of  the  crater  as  a  function  of  depth  during 
the  process  of  formation,  namely 

2  o  ,  2  . 

'  2r^l  +  r^  (22) 

where  r  is  the  crater  radius  at  depth,  1,  and  r^  is  the  radius  of  the  jet.  Secondly, 
final  crater  diameters  in  metals  penetrated  by  shaped  charge  jets  comprised  of  various 
materials  were  compared  with  predictions  of  final  crater  diameter  from  the  theory  for 
final  crater  diameter.  Thirdly,  the  shapes  of  craters  during  formation  in  transparent 
solids  and  liquids  as  photographed  by  means  of  microsecond  framing  camera  techniques 
were  compared  with  predictions  of  the  dynamic  theory. 


V 

■  * 
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Tables  IV,  V  and  VI  present  respectively  comparisons  of  measured  final  crater 
diameters  with  computed  values  for  steel,  copper,  and  aluminum  jets  penetrating 
targets  composed  of  various  metals,  d^  referring  to  the  jet  diameter,  d  the 
computed  final  crater  diameter,  and  d(exp)  the  average  measured  final  crater 
diameter . 


Table  IV:  Comparison  of  computed  and  average  measured  hole  diameters  for  steel  jets 

Jet  striking  velocity  (km/sec) 


Target 

5.1  km/sec(d  =  1.86  mm) 

7  kin/sec(d  =  1.85  tnm) 
o 

8.5  km/sec (d  *=  1.87  mm) 
o 

d' 

d(exp) 

d' 

d(exp) 

d' 

d(exp) 

(mm) 

(mm) 

(mm) 

(mm) 

(mm) 

(mm) 

Lead 

85 

55 

119 

68 

144 

76 

Copper 

20 

21 

27 

27 

34 

32 

Aluminum 

21 

17 

29 

29 

35 

32 

Zinc 

13 

17 

18 

22 

28 

Fe 

13 

13 

18 

16 

22 

21 

Table  V:  Comparison  of  computed  and  average  measured  final  hole  diameters  for 
copper  jets. 

Jet  striking  velocity  (km/sec) _ _ _ 


Target 

5  km/sec (d  =  0.99  mm) 

8.0  km/sec(d^=  0.93  mm) 

8.6  km/sec(d  =  1.61  mm) 
o 

d' 

d(exp) 

d’ 

d(exp) 

d’ 

d(exp) 

(mm) 

(mm) 

(mm) 

(im) 

(mm) 

(mm) 

Lead 

47 

52 

67 

68 

128 

73 

Steel 

7 

9 

11 

15 

20 

20 

Aluminum 

8 

15 

11 

33 

22 

45 

Zinc 

7 

12 

10 

19 

29 

Cu 

11 

11 

16 

16 

31 

31 

Table  VI: 


Target 


Lead 

Copper 

Aluminum 

Zinc 

Steel 


Comparison  of  computed  and  average  measured  final  hole  diameters  for 


aluminum  jets. 

Jet  striking  velocity  (km/sec) 


1  5.4  km/sec(d  =22  mm) 

8.8  km/sec(d  =  1.78  mm) 
o 

9.8  kin/sec(d  =  1.88  mm) 
o 

d' 

d(exp) 

d' 

d(exp) 

d' 

d(exp) 

(mm) 

(mm) 

(mm) 

(mm) 

(mm) 

(mm) 

77 

62 

101 

83 

119 

103 

19 

21 

25 

29 

29 

27 

21 

21 

28 

28 

33 

33 

10 

22 

13 

15 

41 

10 

13 

13 

25 

15 

39 
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In  equation  (3),  which  was  used  to  compute  the  final  crater  diameter,  d^  was 
unknown.  Although  from  flash  radiographs  and  high  speed  photographs  in  vacuum 
obtained  at  this  laboratory  the  jet  diameter  d^  was  known  to  be  of  the  order  of 
magnitude  of  a  millimeter  for  the  particular  shaped  charge  liners  used,  the 
actual  value  of  the  jet  diameter  was  estimated  by  applying  equation  (3)  to  the 
measured  final  hole  diameters  in  targets  of  the  same  metal  as  the  jet  (d^  in  this 
case  being  the  only  unknown).  This  estimated  value  of  the  jet  diameter  was 
subsequently  used  in  equation  (3)  to  compute  the  final  crater  diameters  in  the 
disimilar  metals.  The  comparisons  between  measured  crater  diameter  and  computed 
diameters  were  considered  to  be  reasonably  good  in  view  of  the  fact  there  was 
considerable  uncertainty  in  the  values  of  the  compressive  yield  strengths  involved. 
For  example,  in  the  case  of  lead  values  of  ultimate  tensile  strength  (which  was 
used  in  lieu  of  compressive  yield  strength  since  the  two  are  supposed  to  be 
comparable  and  no  compressive  yield  strength  data  were  available)  have  been  quoted 
which  vary  from  1600  Ib/in"  to  as  high  as  3300  Ib/in  .  In  the  calculations 
2000  Ib/in^  was  used  because  it  was  considered  to  be  most  appropriate  to  the  type 
metal  used  in  the  targets.  However,  had.  3000  Ib/in  been  used  the  agreement  in 
general  would  have  been  much  better.  Table  VII  lists  the  compressive  yield  strengths 
used  in  the  calculations.  Where  compressive  yield  strength  data  were  not  available 
the  ultimate  tensile  strength  was  used  because  the  two  values  are  generally  com¬ 
parable  . 


Material 


Table  VII:  Compressive  yield  or 

(Ib/in^) 


ultimate  tensile  strengths* 


Lead  (cast) 

Steel  (1018,  hot  rolled) 
Copper 

Aluminum  (llOO-H-12) 

Zinc  alloy  (967.  Zn,  47.  Al) 

*From  ASM  Metals  Handbook 


2,000 

70,000 

32,000 

15,000 

40,000  shear 
60,000  Comp 


(dyne /cm  ) 

1.4  X  10® 

5  X  10^ 

2  3  X  10^ 

1  05  X  10^ 
4.2  X  10^ 


% 
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The  formation  of  craters  in  transparent  targets  by  shaped  charge  jets  was 
investigated  by  backlighting  the  target  with  an  explosive  lightbomb  and  photo¬ 
graphing  the  penetration  of  the  Jet  and  the  expansion  of  the  crater  with  a  high 
speed  framing  camera.  The  transparent  substances  studied  were:  water  (  p  =  1.0); 
acetylene  tetrabromide  (  p  =  2.96);  saturated  aqueous  zinc  chloride  solution 
(  p  -  1.88);  gasoline  (  p  *  0.79);  Incite  (  p  -  1.18);  carbon  tetrachloride  (  p  -  1 
and  chloroform  (  p  -  1.49).  The  jet  velocities  used  were  5.1  km/sec,  7.0  km/sec, 
and  8.5  km/sec,  the  jets  in  all  cases  being  steel  and  originating  from  liners 
identical  with  those  used  in  the  final  crater  diameter  studies  in  metals.  Con¬ 
sequently,  the  previously  determined  jet  diameters  were  available  for  use. 

From  the  framing  camera  sequences  the  instantaneous  extent  of  forward 
penetration  was  plotted  against  the  instantaneous  hole  radius  and  the  results 
compared  to  equation  (22).  The  experimental  data  were  found  in  general  to  fall 
quite  well  along  the  theoretical  parabolic  curve.  No  trend  away  from  the 
theoretical  curve  was  discernable  for  gasoline  and  water  or  lucite  and  chloroform 
which  possess  comparable  densities.  However,  in  the  case  of  carbon  tetrachloride, 
zinc  chloride  solution,  and  acetylene  tetrabromide,  each  representing  a  material 
of  increased  density,  there  was  an  increasing  tendency  for  the  crater  diameter  to 
exceed  theoretical  values.  Moreover,  this  tendency  became  more  pronounced  at 
higher  velocities.  At  the  writing  of  AFOSR-TN-60-851 ,  no  explanation  was  advanced 
to  explain  why  the  ratio  of  the  lateral  expansion  velocity  to  the  forward 
penetration  velocity  increased  with  target  density  and  jet  velocity.  However, 
since  that  time  the  thought  has  been  advanced  that  the  effect  could  be  due  to 
increased  compression  of  the  target  material  between  the  crater  and  the  shock 
front  resulting  from  higher  pressure  generated  by  increasing  the  density  of  the 
target  material  or  the  jet  velocity. 

E.  Ionization  and  Electron  Densities  in  Detonation  Solid  Explosives 

The  results  of  this  investigation  were  described  in  item  #1  of  the  list  of 
technical  notes.  This  important  investigation  formed  a  basis  for  some  of  the 
theories  of  detonation  and  external  detonation-generated  plasmas  advanced  later 
by  this  group.  Some  of  the  results  obtained  in  this  study  were  corraborated  by 
papers  given  four  years  later  at  the  Third  ONR  Symposium  on  Detonation,  Princeton, 
New  Jersey.  These  experiments  involved  consisted  of  determining  the  specific 


electrical  conductivities  within  detonation  wavt s  of  solid  explosives  and  mapping 

out  the  shapes  of  the  conducting  regions.  Both  of  these  determinations  were  made 

through  the  application  of  conduction  probe  techniques.  By  means  of  a  theoretical 

expression  for  mobility  of  electrons  under  the  conditions  prevailing  in  detonation 

waves  and  the  measured  conductivities  the  electron  densities  were  calculated.  The 
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resulting  values  which  ranged  from  0.2  x  10  /cm  for  low  density  granular  80/20 

17  3 

AN/TNT  mixtures  to  6  x  10  1/cm  for  Composition  B  were  believed  to  represent 
underestimates  of  the  actual  electron  densities. 

Considerations  were  given  as  to  mechanism  through  which  such  high  electron 
densities  could  be  created.  Thermal  ionization,  for  example,  had  to  be  eliminated 
because  under  the  exceedingly  high  pressure,  relatively  low  temperature  conditions 
involved  in  detonation  waves  the  process  of  thermal  ionization  would  lead  to 
electron  densities  orders  of  magnitude  lower  than  was  measured.  More  or  less 
through  a  process  of  elimination  the  conclusion  was  finally  reached  that  the 
electrons  were  a  result  of  chemionization  which  occurred  during  the  extremely 
rapid  chemical  reactions  involved  in  detonations.  In  the  light  of  this  conclusion 
electron-ion  re- combination  tires  in  the  detonation  wave  were  considered.  Kinetic 
theory  approaches  were  applied  to  the  problem  with  the  result  that  average  re- 
combination  times  were  calculated  to  be  about  10  seconds  for  electron  positive 
ion  recombination  and  they  would,  of  course,  be  shorter  in  the  case  of  negative 
ion  formation.  Such  recombination  times  were  considered  to  be  unrealistically 

short,  however,  because  electron  transition  times  were  known  to  be  the  order 

-8 

of  10  seconds.  On  this  basis  it  was  concluded  that  electron  lifetimes  in  deto¬ 
nating  explosives  were  limited  not  by  purely  kinetic  effects  but  by  recombination 
probabilities  and  that  the  lifetime  was  probably  the  order  of  10  seconds. 

Since  chemical  reaction  times  in  detonation  were  known  to  be  the  order  of 
microseconds  or  roughly  one-hundred  times  the  expected  lifetime  of  an  electron, 
it  was  concluded  that  electrons  would  be  measured  only  in  regions  where  they  were 
being  formed.  Thus,  the  conduction  zone  in  detonation  waves  corresponded  to  the 
chemical,  reaction  zone.  On  the  basis  of  the  high  electron  densities  and  their 
short  lifetimes  it  was  also  concluded  that  at  some  stage  during  the  reaction  each 
molecule  of  explosive  contributes  a  free  electron  and  possibly  an  electron  is 
freed  wherever  a  chemical  bond  was  broken. 
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Measurements  of  the  length  of  the  conduction  zone  were  then  compared  to 
chemical  reaction  zone  lengths  which  were  determined  from  detonation  velocity  vs 

(4) 

diameter  curves  by  the  existing  theories,  namely  the  Eyring  "curved  front"  theory, 
the  Jones  "nozzle"  theory, and  the  "detonation  head"  model  of  Cook.^  ^  These 
comparisons  revealed  that  in  every  case  where  the  explosive  was  ideal  (the  deto¬ 
nation  velocity  equaled  the  "leveled  off"  or  hydrodynamic  value)  the  length  of  the 
conduction  zone  at  the  charge  axis  equaled  within  a  factor  of  two  the  reaction 
zone  length  given  by  the  detonation  head  model,  and  generally  the  agreement  was 
closer  than  this.  Table  VllI  is  a  reproduction  of  a  table  given  in  the  original 
report  showing  these  comparisons.  In  the  case  of  the  two  pon-ideal  explosives, 
coarse  TNT  and  80/20  AN/TNT,  the  conduction  zone  lengths  were  much  shorter  than 
the  reaction  zone  lengths  determined  by  the  "detonation  head"  model. 

The  shapes  of  the  conduction  zones  mapped  out  by  the  probe  technique  were 

then  compared  with  the  "detonation  head"  model.  In  unconfined  cylindrical  charges, 

according  to  the  premises  of  this  model,  a  high  density,  high  pressure  region 

(called  the  detonation  head)  with  approximately  uniform  density  prevails  in  the 

detonation  wave.  When  the  detonation  wave  is  steady  state  and  fully  formed,  this 

high  density  region  is  the  shape  of  a  cone  whose  base  coincides  with  the  detonation 

wave  front  and  whose  sides  are  bounded  by  lateral  rarefaction  or  release  waves 

which  propagate  into  the  charge  starting  at  the  detonation  wave  front.  The  length 

of  the  detonation  head  is  said  to  equal  approximately  its  diameter  which  in  turn 

equals  the  diameter  of  the  charge  minus  a  small  edge  effect  of  about  0.6  cm. 

According  to  this  model  the  detonation  wave  propagates  at  a  velocity  determined 

by  the  extent  of  chemical  reaction  occurring  along  the  streamline  of  maximum 

length  inside  the  detonation  head.  Thus,  explosives  with  a  reaction  zone  length 

a  equal  to  or  less  than  h,  the  length  of  the  detonation  head,  are  ideal,  and 
o 

those  with  longer  reaction  zone  lengths  are  non-ideal.  For  reaction  zone  lengths 
which  were  shorter  than  the  length  of  the  detonation  head,  the  ionized  zones  as 
determined  by  direct  mapping  with  the  probes  were  found  to  have  approximately  the 
shape  of  truncated  cones,  the  length  of  the  zone  becoming  less  toward  the  surface 
of  the  charge.  The  rear  surface  of  the  zone  as  already  stated  corresponded 
roughly  to  the  end  of  the  reaction  zone  predicted  by  the  detonation  head  model  and 
the  sides  of  the  zone  apparently  were  defined  by  the  inward  moving  lateral  rare¬ 
faction  waves  previously  mentioned  that  define  the  detonation  head.  This  is  not 
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surprising  on  the  basis  of  the  assumption  that  the  ionization  is  produced  by 
chemical  reaction  since  it  is  clear  that  the  temperature,  and  consequently  the 
chemical  reaction  rate,  must  drop  sharply  at  the  boundaries  of  the  detonation 
head.  For  non-ideal  explosives  it  was  expected  that  the  conduction  zone  should 
just  correspond  to  the  fully  formed  detonation  head.  In  the  case  of  coarse  TNT 
the  length  at  the  axis  corresponded  well  but  the  off-axis  lengths  tended  to  extend 
past  the  boundaries  of  the  detonation  head.  For  the  80/20  AN/TNT  mixture  the 
conduction  zone  length  at  the  charge  axis  was  somewhat  shorter  than  the  detonation 
head.  This  mixture,  however,  presented  difficulties  which  made  determinations  of 
the  actual  conduction  zone  length  difficult.  With  this  explosive  the  fine  TNT  in 
the  mixture  reacted  much  more  rapidly  than  the  ammonium  nitrate.  After  the  TNT 
reacted,  the  conduction  corresponding  to  the  reaction  of  the  AN  dropped  to  low 
values  which  were. dif ficult  to  measure. 

The  results  of  the  above  study  explained  the  origin  of  the  ionization  in 
external  detonation  generated  plasmas  which  will  be  discussed  later.  They  also 
provided  a  foundation  upon  which  new  theories  of  detonation  were  built.  For 
example, , the  presence  of  such  high  electron  densities  in  detonation  reaction  zones 
with  their  associated  high  mobilities  suggested  that  the  conduction  of  heat  within 
the  wave  was  appreciable  in  contrast  to  the  classical  theory  of  detonation  which  was 
based  upon  transport  phenomena  in  detonation  waves  being  negligible.  In  fact, 
there  are  reasons  to  believe  that  sufficient  heat  conduction  is  a  necessary  con¬ 
dition  for  detonations  to  occur  as  will  be  discussed  later. 

F .  Electrical  Fields  and  Electromagnetic  Radiation  from  Chemical  Detonations 

The  results  of  this  study  were  described  in  item  #7  of  the  list  of  technical 
notes.  During  the  course  of  some  field  work  carried  out  by  this  group,  which 
utilized  the  initiation  of  ejqilosives  by  a  portable  1-1/2  volt  battery  firing  unit, 
the  two  firing  leads  were  accidentally  touched  during  an  explosion,  and  a  strong 
electrical  shock  was  experienced.  A  further  study  of  this  phenomenon  revealed  that 
appreciable  electrical  fields  existed  around  explosions.  These  preliminary  studies 
provided  motivation  to  investigate  the  types  of  electrical  manifestations  resulting 
from  chemical  detonations  and  the  mechanisms  whereby  they  were  generated. 

Due  to  the  fact  that  such  a  large  number  of  varied  experiments  were  carried  out 
as  a  part  of  the  program  it  is  not  practical  to  describe  them  all  in  this  final 
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report.  Hence,  only  the  main  general  conclusions  will  be  discussed  as  well  as  a 
brief  description  of  the  general  experimental  method. 

Measurements  were  made  for  charges  ranging  in  size  from  lOg  to  approximately 
25  kg  over  a  frequency  range  extending  from  a  few  cycles/sec  to  500  megacycles. 
Most  of  these  measurements  were  made  with  the  charge  to  pickup  separation  ranging 
from  several  meters  to  a  hundred  meters  or  so,  which  would  lie  well  within  the 
induction  zone  for  frequencies  below  the  megacycle  range.  It  was  expected  at 
first  that  most  of  the  signals  would  be  above  10  kc.  Hence,  in  first  experiments 
a  broad  band  system  consisting  of  a  2  meter  length  whip  antenna  feeding  directly 
into  the  pre-amplifier  stage  of  a  Tektronix  545  oscilloscope  was  used.  Results 
soon  showed,  however,  that  the  strongest  signals  were  in  a  lower  frequency  range, 
which  dictated  that  equipment  covering  a  broader  spectrum  be  used.  For  subsequent 
measurements,  depending  on  the  information  desired,  one  or  more  of  the  following 
pickups  were  used:  a  5000  megohm  input  Impedance  cathode  follower  with  a  maximum 
frequency  response  of  100  kc,  a  1000  megohm  input  Inpedance  cathode  follower,  a 
1  megohm  input  impedance  completely  transistorized  emitter-follower  with  a  high 
frequency  cut  off  of  1  megacycle,  a  500  me  receiver,  a  98  cm  receiver,  and  a  50  me 
receiver.  The  cathode  followers  were  all  battery  operated  and  were  designed  such 
that  their  pickup  antennae  were  maintained  at  ground  potential,  10  cm  diameter 
copper  spheres  being  used  for  antennae.  In  each  case  the  low  frequency  cut  off 
was  determined  by  the  input  impedance,  the  distance  from  the  signal  source,  and 
the  type  antenna  employed.  A  horizontal  yagi  antenna  was  used  with  the  500  me 
receiver,  while  with  the  98  me  and  50  me  receivers  one-half  wave  length  vertical 
or  horizontal  dipole  antennae  were  used.  In  all  cases  recordings  were  made  with 
Tektronix  545  oscilloscopes. 

A  10  KV  square  wave  calibrator  was  constructed  to  assess  the  linearity  of 
response  of  the  low  frequency  cathode  follower  systems  at  frequencies  essentially 
down  to  10  cps.  This  cdlbrator  charged  a  25  cm  diameter  spherical  conductor  with 
an  8  KV  square  wave,  the  signal  from  which  was  picked  up  by  the  cathode  follower 
systems.  It  was  found  possible  to  achieve  nearly  linear  D.C.  response  for  6m 
separation  between  the  pickup  and  the  25  cm  ball  by  use  of  a  0.01  mfd  capacitor 
connected  from  the  antenna  to  ground  in  the  case  of  the  1000  megohm  input 
impedance  cathode  follower  system.  Thus,  when  very  low  frequency  response  was 
desired  such  an  integrating  capacitor  was  used.  A  capacitance  of  0.01  mfd  was  a 
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compromise  value,  however,  as  larger  capacitors  lead  to  somewhat  more  linear 
response  but  they  decreased  the  gain  prohibitively.  The  rise  time  of  the  1000  meg¬ 
ohm  input  impedance  cathode  follower  system  under  these  circumstances  was  found  to 
be  about  10  ijsec  from  which  an  estimate  of  50  KC  for  the  high  frequency  cut  off  was 
obtained . 

One  of  the  greatest  difficulties  encountered  in  carrying  out  the  experiments 
was  the  reduction  of  background  noise  to  useable  levels.  This  noise  resulted 
primarily  from  circuitry  employed  originally  to  trigger  the  oscilloscope  and  to 
fire  the  electric  detonators  which  initiated  the  charges  The  first  problem  was 
circumvented  by  employing  a  battery  powered  transistorized  photocell  trigger  unit 
which  triggered  the  oscilloscopes  from  the  first  light  either  of  the  explosion  or 
an  auxiliary  flash  which  was  timed  to  occur  a  few  microseconds  before  the  explosion. 
Electric  detonator  systems  were  abandoned  early  in  the  program  because  the  presence 
of  wires  near  the  explosion  was  found  to  alter  the  signals  considerably.  Experience 
also  proved  that  any  obstructions  rear  the  charge  exerted  an  influence  on  the 
signals.  Hence,  in  most  tests  the  charges  were  suspended  above  the  ground,  away  from 
any  obstructions,  by  a  nylon  string  In  order  that  the  charge  be  as  free  as  possible 
from  external  Influences  at  the  time  of  detonation  a  special  method  of  initiation 
was  employed.  An  ordinary  fuse  blasting  cap  was  used  to  initiate  the  charge,  the 
fuse  being  ignited  electrically  by  a  squib.  As  near  to  the  blasting  cap  ac  possible 
(about  1")  the  fuse  was  cut  and  a  section  of  JPN  propellant  was  inserted  which  upon 
"burning  through"  allowed  the  burned  fuse  to  separate  and  the  section  connected  to 
the  squib  and  the  squib  lead  wires  to  fall  to  the  ground.  Thus,  before  the  explosion 
occurred  the  charge  hung  perfectly  free. 

The  more  important  experimental  results  that  were  obtained  may  be  summarized 
briefly  as  follows: 

(1)  The  largest  voltages  generated  by  chemical  detonations  were  found  to  lie 
in  a  frequency  range  below  a  few  hundred  cycles /second .  Only  a  small  fraction  of 
the  low  frequency  signals  was  radiated. 

(2)  The  dominant  frequency,  provided  the  explosive  was  free  from  obstructions, 
was  found  to  decrease  with  charge  size.  Also,  the  electrical  potential  generated 

by  an  explosion  was  found  to  increase  with  charge  size.  Accurate  absolute  evaluations 
appeared  to  be  impossible  because  the  geometry  and  charge  distribution  of  the  signal 
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source  created  by  the  explosion  was  not  accurately  known.  However,  in  magnitude  the 
potential  appeared  to  Increase  in  direct  proportion  to  charge  weight.  By  means  of 
direct  calibrations  of  equipment  and  assuming  that  the  signal  source  was  equivalent 
to  an  expanding  charged  sphere  whose  diameter  equaled  the  diameter  of  the  deto¬ 
nation  products,  -18,000  V  was  determined  to  be  the  maximum  signal  source  potential 
for  a  552g  pentolite  charge. 

(3)  Accurate  determinations  of  attenuation  with  distance  of  the  detonation 
generated  electrical  fields  were  difficult  because  the  frequency  response  of  the 
equipment  varied  with  the  charge  to  pickup  distance. 

(4)  As  a  rule  initially  the  electrical  potential  of  the  detonation-generated 
signal  source  became  uniformly  more  negative  with  time  until  a  maximum  was  reached. 
Thus,  as  charge  was  lost,  its  potential  became  less  negative.  In  most  instances 
the  signal  source  potential  remained  negative  with  respect  to  the  earth  for  the 
duration  of  the  signal.  Intermediate  frequencies  seemed  to  be  generated  largely 
by  the  expanding  cloud  of  detonation  products  hitting  obstructions  or  contacting 
grounded  conductors 

(5)  Ihe  reproducibility  of  the  detonation-generated  electrical  fields  was 
markedly  affected  by  obstructions,  principally  grounded  conductors  (even  though 
they  possessed  high  Impedances)  which  were  contacted  by  the  expanding  detonation 
products.  Another  even  more  important  factor  was  the  blowing  of  dust  from  the  earth's 
surface  by  the  explosion.  Enclosing  the  charge  by  a  layer  of  plaster  of  paris, 

for  example,  markedly  increased  the  signals  at  all  frequencies.  The  low  frequency 
component  could  be  enhanced  by  a  factor  as  high  as  20  by  this  process. 

(6)  Increased  outputs  from  cathode  followers  covering  both  a  low  frequency 
range  and  an  intermediate  frequency  range  were  observed  when  the  charge  height 
was  increased  within  certain  limits.  Such  increases  did  not  result  from  an 
Increased  potential  of  the  detonation-generated  signal  source  but  rather  from 
decreased  distortion  of  the  electrical  field  by  the  earth  as  the  height  was 
increased.  This  fact  was  borne  out  by  static  calibrations  in  the  laboratory. 

(7)  No  significant  effects  of  charge  shape,  charge  orientation  or  position  of 
the  detonator  on  the  charge  were  observed  on  low  frequency  and  intermediate  frequency 
electrical  fields.  The  effect  upon  the  high  frequencies  in  this  regard  was  not 
determined . 
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(8)  Experiments  were  carried  out  in  a  RF  cage  which  was  sufficiently  large  to 
contain  both  the  explosion  and  the  pickup  antennae.  Means  were  available  to  control 
over  a  wide  range  the  static  electrical  field  within  the  cage.  Results  showed  that 
the  static  electrical  field  surrounding  an  explosion  had  little  effect  upon 
explosion-generated  electrical  fields  or  the  EM  radiation. 

(9)  The  high  frequency  radiated  energy  occurred  in  the  form  of  short,  random 
bursts  which  covered  a  broad  band  of  frequencies.  In  many  cases  some  of  the  bursts 
of  radiation  were  observed  to  correspond  to  higher  frequency  transients  manifested 
on  oscillograph  traces  of  the  outputs  from  low  frequency  response  cathode  followers. 

(10)  The  detonation- generated  radiation  bore  a  strong  resemblance  to  radiation 
resulting  from  corona  discharge. 

(11)  With  bare  charges  there  was  found  to  be  a  non-reproducible  delay  between 
initiation  and  the  first  bursts  of  radiation.  This  delay  was  generally  at  least 
several  hundred  microseconds  and  could  be  as  much  as  several  milliseconds.  When 
charges  were  enclosed  by  a  layer  of  plaster  of  paris  this  delay  was  reduced. 

(12)  Increasing  the  charge  size  increased  primarily  the  number  of  bursts  of 

radiation  rather  than  their  amplitude. 

(13)  Most  of  the  radiation  pulses  contained  both  vertical  and  horizontal 
electric  field  components,  although  cases  were  observed  where  the  bursts  were 
completely  vertically  or  horizontally  polarized. 

(14)  During  the  later  stages  of  the  expansion  process,  when  the  explosion  was 
sufficiently  close  to  the  earth,  the  vertical  component  of  the  radiation  was  found 
to  predominate  and  was  markedly  enhanced  at  the  instant  of  impact  of  the  gas  cloud 
with  the  ground. 

On  the  basis  of  the  above  results  a  mechanism  was  proposed  as  follows:  The 
results  cited  concerning  the  large  enhancement  of  signal  when  the  charges  were 
enclosed  by  inert  solids  show  that  the  air  shock  generated  by  the  blast  or  the 
external  explosive  plasma  are  unimportant  as  sources  of  electrical  signals  because 
a  layer  of  the  inert  material  eliminates  the  plasma  and  attenuates  the  air  shock 
considerably.  The  enhancement  of  the  low  frequency  signals  produced  by  inert, 
solid  confinement  can  be  visualized  to  occur  in  the  following  way:  The  force  of 
the  explosion  accelerates  to  high  velocities  the  small  particles  formed  when  the 
solid  is  pulverized.  These  small,  high  velocity  particles  pick  up  an  electrokinetic 


charge  as  they  move  through  the  atmosphere  thus  forming  an  expanding  cloud  of 
charged  particles.  By  virtue  of  the  fact  that  these  inert  particles  pick  up  more 
and  more  charge  as  time  goes  on  and  at  the  same  time  the  expanding  cloud  grows 
larger  and  larger,  the  electrical  field  sensed  by  the  pickup  increases.  The  length 
of  time  that  the  electrical  field  continues  to  grow,  as  well  as  its  magnitude,  would 
therefore  be  expected  to  increase  with  charge  size.  This  external  field  finally 
decays  to  zero  when  the  velocity  of  the  J>ar tides  become  so  low  that  they  lose  their 
charge.  Contact  of  the  expanding  cloud  with  the  ground  will  Increase  the  rate  at 
which  it  loses  its  charge  as  shown  by  charge  height  studies.  Contact  with  the 
ground,  other  conductors  and  probably  also  even  other  rather  massive  nonconductors 
can  introduce  transients  in  this  external  field.  These  transients,  however,  may 
or  may  not  be  abrupt  enough  to  generate  megacycle  components  which  may  be  measured 
as  radiation.  In  fact,  the  electrokinetic  effect  of  small,  high  velocity,  solid 
particles  may  also  account  for  the  observed  electrical  effects  even  for  bare  charges. 
All  explosives  exhibit  an  "edge  effect."  That  is,  a  small  layer  of  explosive  around 
the  charge  is  not  immediately  consumed  by  the  detonation  in  unconfined  charges. 
Possibly  some  of  this  material  may  be  propelled  outward  in  the  form  of  high  velocity 
microparticles.  However,  this  explanation  for  bare  charges  is  not  justified  from 
the  evidence  yet  available 

As  mentioned  above,  transients  in  the  external  field  can  be  produced  in  various 
ways.  It  is  felt  that  the  high  frequency  components  (i.e.,  50  MC,  98  MC  and  500  MC) 
may  be  generated  primarily  by  miniature  "lightning  bolts"  occurring  between  regions 
of  differing  potential  within  the  signal  source.  This  suggestion  seems  to  be  borne 
out  by  the  observed  delays  between  the  detonation  and  the  first  bursts  of  radiation 
as  well  as  the  randomness  of  the  bursts  and  their  resemblance  to  radiation  generated 
by  corona  discharge.  This  model  to  account  for  high  frequencies  also  seems  to 
agree  with  the  effect  of  charge  size.  That  is,  there  appears  under  a  given  set  of 
conditions,  to  be  a  minimum  charge  size  which  will  generate  frequencies  in  the  mega¬ 
cycle  range.  As  charge  size  is  increased,  the  most  striking  effect  produced  is 
that  the  radiation  pulses  occur  much  more  frequently  with  only  a  small  increase  in 
magnitude.  This  suggests  that  a  critical  value  of  potential  gradient  must  exist 
within  the  expanding  cloud  in  order  for  a  pulse  to  be  generated. 
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It  has  been  noted  that  when  charge  sizesaround  lOOg  are  used  no  50  MC  or  98  MG 
signals  are  detected,  while  at  other  times  using  the  same  sensitivity  settings  signals 
will  be  received.  This  suggests  that  some  properties  of  the  surrounding  gaseous 
medium  subject  to  change  may  determine  whether  or  not  critical  conditions  may  attain 
within  the  signal  Source  cloud  One  possible  factor  is  humidity.  However,  no 
attempts  have  so  far  been  made  to  evaluate  this  factor 

G.  External  Detonation-Generated  Plasmas 

The  results  regarding  this  subject  have  been  described  in  items  #5,  #10  and 
#12  of  the  list  of  technical  notes 

When  the  detonation  wave  front  reaches  a  free  surface  of  a  charge  in  air,  for 
example  the  end  of  a  cylindrical  charge,  a  highly  luminous  cloud  is  formed  which 
propagates  away  from  the  charge  at  high  velocity.  This  phenomenon  had  been  observed 
by  others  as  well  as  members  of  this  group  for  a  long  time  but  had  attracted  little 
Interest  because  the  luminosity  was  thought  to  be  simply  a  result  of  thermal 
ionization  In  a  shock  wave  generated  in  the  atmosphere  surrounding  the  charge. 

While  studying  microsecond  framing  camera  photographs  in  color  showing  detonations 
of  liquid  explosives,  it  was  noted  that  the  luminous  cloud  which  appeared  at  the 
free  surface  of  the  liquid  exhibited  a  remarkable  tendency  to  "hang  together"  as  it 
propagated  through  the  atmosphere.  The  Interest  kindled  by  this  observation  led  to 
further  investigations  of  the  phenomenon  which  were  carried  out  under  this  contract. 

Through  probe  conduction  techniques,  these  luminous  clouds,  as  expected,  were 

found  to  be  highly  ionized.  The  electron  densities,  in  fact,  were  estimated  to  be 
16  3 

the  order  of  10  1/cra  or  about  1/10  to  1/100  the  density  prevailing  in  detonation 
reaction  zones.  The  first  concern  naturally  was  determining  the  origin  of  this 
ionization  Thus,  the  mechanisms  which  could  possibly  be  effective  were  formulated 
and  considered.  These  mechanisms  were' 

(1)  Unreacted  or  partially  reacted  explosive  is  spalled  from  the  end  of  the 
charge  upon  the  arrival  of  the  detonation  wave.  This  material  reacts  as  it  is 
propelled  from  the  charge  and  by  virtue  of  this  chemical  reaction  forms  ionized 
species . 

(2)  The  luminous  cloud  consists  largely  of  ions  and  electrons  propelled  from 
the  detonation  reaction  zone.  As  a  result  of  chemical  reaction  occurring  at 

the  surface  of  the  explosive  when  the  detonation  wave  emerges,  ions  and  electrons 
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»hlch  are  fomed  are  propelled  forward  Into  a  ™eh  lower  den.lty  region,  the 

ejected  Ufettae  of  an  electron  therefore  being  nntch  longer  hy  virtue  of  the  reduced 
density  of  gas  molecules. 

(3)  Ihemal  lonltatlon  within  a  ehock  wave  that  la  generated  In  the  gaseous 
medium  outside  the  charge  upon  the  emergence  of  the  detonation  wave 

Mechanisms  (1)  and  (3)  were  soon  discounted  hy  early  experiments,  and  this 

conclusion  was  later  verified  hy  a  large  number  of  facts  gained  through  subsequent 

research^  garly  experiments,  for  example,  were  performed  In  which  the  lonltaLon 

generate  by  pure  shock  waves  was  compared  to  that  produced  when  a  detonation  wave 

g  rom  the  end  of  a  bare  charge.  Among  these  comparisons  possibly  the  most 

sinking  was  that  obtained  using  a  charge  comprised  of  a  92/8  mixture  of  fine 

explosive  grade  ammonium  nltr.te-nltropropane  whose  end  was  bare  and  a  2"dla  x  8" 

length  Composition  B  charge  whose  end  was  terminated  with  a  1/V  thick  glass  plate. 

8  ass  ncldentally  being  a  fairly  good  Impedance  match  tor  Composition  B.  In  the 

case  of  the  glass-termlnated  Composition  B  charge  no  luminosity  was  apparent  on  a 

ram  ng  camera  sequence  of  the  event.  However,  a  luminous  cloud  propagated  from 

the  tree-ended  anmx.nlum  nltr.te-nltropropane  charge  (which  was  photographed 

slmu  taneously  with  the  Composition  B  detonation)  when  the  detonation  wav.  reached 

e  re.  end  Calculations  revealed  that  the  pressure  of  the  shock  wave  driven  by 

the  glass-termlnated  Composition  B  charge  was  at  least  as  great  as  the  pressure  of 

a  shock  wave  which  could  be  driven  by  the  b.re-ended  AN-HP  charge.  Such  detonation- 

generated  luminous  clouds.  In  tact,  have  been  observed  to  propagate  from  tree 

sur  aces  of  charges  no  ,»tter  how  weak  the  detonation  Is  as  long  as  It  occurs 

Respectable  luminous  clouds,  for  example,  have  been  observed  to  be  generated  by 

etonatlons  In  99/6  ammonium  nitrate-fuel  oil  mixtures  where  the  detonation  pressure 
was  as  low  as  13  kb.  ^ 

As  mentioned  earlier,  framing  camera  sequences  obtained  of  detonation  liquid 
explosives  In  glass  beakers  showed  that  the  luminous  cloud  generated  at  the  surface 
of  the  liquid  did  not  expand  significantly  as  It  left  the  top  of  the  beaker,  which 
is  not  the  behavior  one  would  expect  from  a  shock  wave.  Furthermore,  the  cloud 
did  not  exert  sufficient  external  pressure  to  fracture  the  beaker  before  the  deto¬ 
nation  products  following  somewhat  behind  shattered  It.  Experiments  with  conduction 
probes  revealed  an  Interesting  pulsating  phenomenon.  Under  certain  conditions  the 
length  of  the  luminous  cloud  was  observed  to  oscillate  and  the  conduction  varied  with 
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the  oscillations,  being  the  greatest  when  the  plasma  length  was  the  shortest  and 
the  least  when  it  was  the  longest. 

On  the  basis  primarily  of  the  aforementioned  evidence,  mechanism  (2)  was 
considered  the  most  app licable,  and  the  luminous  ionized  clouds  were  termed  "external 
detonation-generated  plasmas"  in  contrast  to  the  ionization  in  the  detonation  reaction 
zone  which  was  termed  the  "internal  detonation-generated  plasma".  .Since  such  high 
electron  and  ion  densities  were  known  to  exist  in  detonation  reaction  zones,  it  was 
not  surprising  that  some  evidence  of  this  would  be  manifested  exterior  to  the  charge 
when  the  detonation  wave  emerged.  (Considerable  evidence  to  support  this  conclusion 
was  obtained  during  subsequent  studies,  the  more  important  points  will  be  given 
below  )  In  order  to  account  for  the  remarking  cohesive  properties  of  the  external 
detonation-generated  plasmas  the  "quasi-lattice"  or  "metallic  model"  for  the  plasma 
was  advanced  According  to  this  model  the  electrons  circulate  through  the  positive 
ion  lattice  in  correlated  patterns  similarly  as  occurs  in  metals.  In  other  words, 
the  external  detona tion- generated  plasmas  were  considered  to  involve  the  formation 
of  the  outermost  valence  electrons  from  localized  states  in  atoms  or  molecules  to 
collective  electron  states  in  which  they  moved  freely,  but  at  an  appreciable  negative 
potential  energy  throughout  the  plasma  to  give  it  appreciable  cohesion  The 
existence  of  a  quasi-lattice  state  in  external  detonation-generated  plasmas  inferred 
the  existence  of  a  similar  state  for  the  internal  detonation-generated  plasmas,  or, 
in  other  words,  that  the  detonation  state  was  a  "plasma"  state.  Such  a  structure 
would  explain  the  exceedingly  high  electron  densities  in  detonation  reaction  zones 
because  the  electrons  would  not  be  completely  free  but  would  still  be  bound  by  the 
cohesive  energy  of  the  lattice,  and  therefore  the  state  would  be  easier  to  achieve 
than  a  completely  free  electron  state  According  to  the  quasi - lattice  model  of  the 
plasmas,  the  greater  the  electron  density,  the  greater  is  the  negative  potential 
energy  of  the  system  and  correspondingly  the  greater  is  the  cohesion  of  the  plasma. 
However,  the  quasi- latt Ice  is  opposed  by  thermal  energy  As  long  as  the  conditions 
prevail  that  E  »  kT  (E^  being  the  cohesive  energy),  then  the  plasma  will  be 
stable,  but  when  E^  decreases  to  approach  kl  then  the  plasma  will  decay  to  an 
ordinary  gaseous  state 

On  the  basis  of  mechanism  (2)  the  initial  length  of  the  external  detonation¬ 
generated  plasma,  before  the  plasma  has  time  to  expand  significantly,  should  be 
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proportional  to  the  reaction  time  of  the  explosive,  i.e.,  L*  -  kV^-r  where  L*  is 
the  initial  plasma  length,  k  is  a  proportionality  constant,  V*  is  the  Initial  plasma 
velocity,  and  t  is  the  reaction  time.  Measurements  made  with  nitromethane,  Dithe- 
kite  13  and  tetryl  revealed  that  such  was  indeed  the  case  with  k  equaling  approximately 
0.5.  Further,  measurements  with  loose  tetryl  of  varying  particle  size  revealed  that 
the  initial  external  detonation-generated  plasma  length  was  directly  proportional  to 

the  particle  size  as  is  the  reaction  time. 

Subsequent  work  vielded  a  great  deal  of  information  concerning  such  things  as 


the'  origin,  structure,  decay  characteristics,  etc.,  of  external  detonation-generated 
plasmas.  The  more  significant  findings  are  summarized  briefly  below: 

(1)  The  decay  of  electrons  in  both  internal  and  external  plasmas  takes  place 
primarily  by  the  formation  of  negative  ions  with  subsequent  neutralization  of  the 
negative  ions  and  positive  ions.  However,  some  direct  positive  ion  capture  of 
electrons  would  be  expected  to  occur.  These  decay  reactions  are  represented  by 
the  following  equations 


(not  necessarily  equilibrium) 
+€  -T  b”  (not  necessarily  equilibrium) 
B*  +  -»  neutral  products 


(23) 

(24) 

(25) 


The  importance  of  negative  ion  formation  in  the  decay  of  plasmas  was  illustrated 
by  decay  measurements  performed  for  plasmas  propagated  into  various  gaseous 
atmospheres,  the  most  rapid  decay  being  observed  for  gases  which  possessed  the 
highest  electron  affinities  This  result  showed  also  that  the  gas  through  which 
the  plasma  was  propagated  was  an  important  contributor  to  the  reduction  of  the 
electron  density  in  the  plasma  and  its  subsequent  decay  to  a  gas.  The  electrical 
conductivities,  measured  for  plasmas  propagated  into  various  atmospheres,  were 
extrapolated  to  the  charge  surface  with  the  result  that  the  initial  conductivities 
possessed  the  same  values,  indicating  again  that  the  major  share  of  the  ionization 
originated  from  the  explosive  rather  than  being  generated  in  the  gas  surrounding 

the  charge 

(2)  The  conduction  zone  of  the  external  plasma  corresponds  precisely  to  its 
luminosity  zone.  This  supports  the  concept  that  the  luminosity  is  a  result  of 
quantum  transitions  and  conduction  is  proportional  to  the  concentration  of  free 


electrons . 
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(3)  External  detonation- generated  plasmas  are  actually  relatively  high  in 
density  and  therefore  propagate  with  a  relatively  large  amount  of  material  in  them. 
This  fact  was  illustrated  by  framing  camera  pictures  showing  the  decay  and 
"explosions"  of  plasmas  to  a  gas.  The  volume  of  gas  produced  as  a  result  of  such 
explosions  has  been  estimated  to  be  as  much  as  100  times  the  original  plasma  volume. 
The  decay  from  the  plasma  state  to  the  gaseous  state  apparently  occurs  when  the 
cohesive  energy  of  the  plasma  as  a  result  of  loss  of  electrons  drops  to  a  value 
approximating  the  thermal  energy.  One  would  expect  that  plasmas  with  high  electron 
densities  such  as  produced  by  fast  reacting,  high  density,  high  velocity  explosives 
would  be  more  stable  than  those  from  "slower  explosives,  and  this  has  proven  to  be 
the  case.  The  decay  or  explosion  of  plasmas  to  a  gas  is  promoted  by  interactions 
which  tend  to  "stir"  the  plasma  and  weaken  its  lattice  structure  such  as  rapidly 
compressing  a  plasma  and  then  allowing  it  to  expand  or  impacting  a  barrier. 

Under  proper  lighting  conditions  a  non-luminous  shock  wave  may  be  observed  to 
separate  from  the  plasma  when  it  explodes.  The  interpretation  of  this  phenomenon 
is  that  when  the  electron  density  in  the  plasma  has  dropped  sufficiently  that  the 
plasma  state  decays  to  the  gaseous  state,  the  ensuing  rapid  expansion  of  the  gas 
drives  a  shock  ahead  of  it. 

(4)  External  plasmas  may  decay  in  one  gas  and  later  reappear  in  another, 

e.g.,  they  have  been  shown  to  disappear  in  propane  and  reappear  in  air.  This 
result  may  be  explained  on  the  basis  that  the  ion  decay  is  small  via  equation  (25) 
but  the  reactions  of  the  type  +  B^  are  fast.  Now  B^  ions  from  propane 

do  not  ionize  rapidly  and  therefore  serve  as  traps  for  electrons.  However,  the 
do  ionize  rapidly.  Therefore,  as  soon  as  the  propane  medium  becomes  diluted  with 
air  the  free  electrons  reappear. 

(5)  Plasmas  are  readily  manipulated  by  mechanical  means.  They  can  be  com¬ 
pressed  from  large  diameters  into  much  smaller  diameter  glass  tubes,  They  can  also 
be  readily  "piped"  around  corners.  In  such  experiments  the  cracking  of  the  glass 
tubing  was  not  sufficient  to  prevent  transparency  until  a  microsecond  or  so  after 
the  plasma  had  gone  by. 

(6)  Manipulations  with  magnetic  fields,  though  attempted,  have  not  as  yet 
proved  successful  presumably  because  of  insufficient  magnetic  field  intensity. 


( 
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(7)  Framing  camera  sequences  of  softly  backlighted  plasmas  from  Dithekite-13, 
which  is  oxygen  balanced  and  possesses  nearly  transparent  detonation  products,  have 
shown  some  very  interesting  structural  features  of  external  detonation-generated 
plasmas;  a  diagram  of  one  picture  from  a  sequence  is  given  in  Fig.  2.  The  sequence 
clearly  showed  a  complicated  multiple  shock  pattern  which  one  would  not  conceive 
of  existing  were  the  plasma  merely  shock  induced  ionization,  because  under  this 
circumstance  only  a  strong  shock  at  the  plasma  front  would  exist.  The  resulting 
shock  pattern  generated  in  the  Incite  witness  block  by  impact  of  the  plasma  showed, 
however,  that  the  shock  corresponding  to  the  front  of  the  plasma  was  weak  in  comparison 
to  the  shock  at  the  rear  of  the  plasma.  As  already  described  during  their  existence 
of  several  hundred  microseconds  or  so,  if  undisturbed  in  flight,  external  detonation¬ 
generated  plasmas  undergo  accelerating  decay  which  sometimes  reaches  explosive 
magnitudes  Decay  occurs  most  rapidly  in  regions  of  low  electron  density  and  the 
decay  promotes  further  decay.  Since  the  electron  densities  are  highest  at  the 
front  of  the  plasma  and  lowest  at  the  rear  (as  determined  by  probe  measurements) 
the  decay  occurs  most  rapidly  at  the  rear.  Even  though  the  density  of  a  plasma 


backlighted 
Incite  block 

jet  from  plasma 


plasma  bow  wave 

recombination 
shock  of  plasma 

bow  wave  corre¬ 
sponding  to 
recombination  shock 

Figure  2  Diagram  of  one  picture  from  a  framing  camera  sequence  of  a 
"softly"  backlighted  plasma  from  dltheklte  13. 
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may  be  several  hundred  times  that  of  air  the  external  detonation-generated  plasmas 
exhibit  little  external  pressure  owing  to  their  high  cohesion  When,  however,  decay 
of  the  plasma  to  a  gas  occurs,  the  pressure  rises  sharply.  Therefore,  at  the  rear 
of  the  plasma  where  decay  is  occurring  rapidly,  a  rapid  pressure  rise  should  occur. 

The  decay  and  sudden  rise  of  pressure  at  the  rear  of  the  plasma  is  believed  to  be 
the  cause  of  the  strong  shock  at  the  rear  of  the  plasma  For  this  reason  this 
shock  was  designated  the  re-combination  shock.  The  weak  shock  at  the  plasma  front, 
on  the  other  hand,  was  considered  to  be  simply  a  compressional  bow  wave  produced  by 
the,  plasma  propagating  through  the  atmosphere  According  to  this  model  external 
detonation-generated  plasmas  would  be  "jet  propelled"  by  decay  of  the  plasma  at 
the  rear  to  a  gas,  Indeed,  the  jetting  of  gases  off  the  back  of  the  plasma  was 
suggested  in  the  framing  camera  seque<^ce  from  which  Fig.  2  was  taken,  what  appeared 
to  be  a  Mach  diamond  being  clearly  evident,  It  is  perhaps  this  jet  propulsion  in 
external  detonation-generated  plasmas  that  causes  them  to  separate  from  the  main 
detonation  products  and  propagate  generally  at  super  detonation  speeds.  Admittedly, 
the  above  explanation  departs  widely  from  classical  theory  However,  there  is  no 
way  evident  to  us  whereby  the  classical  theory  of  thermal  ionization  in  a  shock 
wave  can  explain  the  complicated  shock  patterns  shown  in  Fig  2.  These  patterns 
we  feel  are  a  justification  for  the  model  of  the  plasma  which  has  been  proposed, 

(8)  When  external  detonation-generated  plasmas  were  conpressed  the  conductivity 
was  found  to  increase  in  a  manner  comparable  to  that  in  metals  under  compression. 

(9)  A  more  extensive  study  involving  plasmas  propagated  into  a  wide  variety 
of  gases  was  carried  out,  the  list  of  which  included  argon,  air,  nitrogen,  oxygen, 
carbon  monoxide,  carbon  dioxide,  ammonia,  chlorine,  propane,  helium,  and  hydrogen, 
all  under  atmospheric  pressure  and  in  addition  hydrogen  and  air  at  about  0.5  mm 
pressure.  These  studies  in  agreement  with  the  less  extensive  earlier  work  involving 
plasmas  propagated  into  a  number  of  gases  showed  that  magnitude  of  the  conductivity 
extrapolated  to  the  charge  surface  was  essentially  independent  of  the  gas, and  the 
decay  rate  of  the  plasma  was  strongly  dependent  upon  the  gas  as  was  the  velocity  of 
the  plasma.  In  addition,  it  was  found  that  the  plasma  could  be  luminous,  transparent, 
or  dark  depending  upon  the  gaseous  medium  and  its  pressure,  as  shown  by  the  results 
given  in  Table  IX.  Measurements  of  conductivity  therefore  were  concluded  to  be  the 
only  reliable  means  of  detennining  whether  or  not  a  plasma  existed.  It  was  interesting 
to  note  that  even  in  the  case  of  a  vacuum  ionized  clouds  (plasmas)  were  propelled  from 
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Table  IX:  Qualitative  comparisons  of  visible  light  intensity  for  ionization  waves 
propagated  through  various  gases  at  ambient  pressure. 


Explosive 

Gaseous  atmosphere 

Ionization 

potential 

(e.v.) 

Visible  luminosity 

D-13* 

Argon 

15.7 

Extremely  brilliant 

D-13 

Air 

-- 

Bright 

Nitrogen 

15.5 

Bright 

Oxygen 

12.5 

Bright 

Carbon  monoxide 

14.1 

Bright 

Carbon  dioxide 

14.5 

Bright 

D-13 

Aitmonia 

11.2 

Moderately  bright 

Chlorine 

13.2 

Moderately  bright 

D-13 

Propane 

-  - 

Nonluminous-dark  cloud 
forms 

D-13 

He  1 ium 

24.5 

No  detectable  luminosity 
transparent  cloud  only 

Hydrogen 

15.6 

Same  as  in  He 

Composition  B  and  tetryl 

He 1 ium 

24.5 

Weak  momentary  light  at 
charge  surface.  No 
luminosity  at  long 
distances  from  charge 

Composition  B 

Air  vacuum 

Very  low  to  no  luminosity 
depending  on  actual 
pressure 

*  D-13  =  63/24/13  HNO^/ni trobenzene/water 

free  surfaces  of  detonating  explosives.  Experimental  measurements  were  compared  to 
calculations  of  electron  density  based  upon  hydrodynamic  theory  and  the  mechanism 
of  thermal  ionization  in  shock  waves  which  were  assumed  to  be  traveling  at  the 
experimentally  determined  plasma  velocities.  These  comparisons  were  carried  out 
for  both  air  and  hydrogen  at  initial  pressures  of  645  mm  Hg  and  0.5  mm  Hg.  In  order 
to  minimize  ambiguities  that  might  result  from  the  use  of  an  incorrect  probe 
calibration  factor  required  to  transform  conductance  data  to  conductivity  data, 
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ratios  of  the  conductivities  at  645  mm  Hg  and  0.5  mm  Hg  were  also  compared.  In  the 
case  of  air  near  the  surface  of  a  bare  Composition  B  charge  o*  -  0.012  mho/cm 

and  0  Q  3  =  0.006  mho  cm^were  measured  while  -  29  and  -  130  mho/cm  were 

calculated  giving  o  545/0^^5  -  2  and  3  =  0.2.  When  the  plasma  was  extruded 

into  a  3/8"  diameter  plastic  tube,  the  results  were  =  0.52  mho/cm,  o*^  3  -  0.26 

mho/cm,  o^^3  =  130  mho/cm  and  3  =>  160  mho/cm  giving  the  ratios  0*545/0*0  5  =  2"0 
°645^°0.5  “  Thus,  for  air  the  disagreement  was  not  sufficient  to  be  con¬ 

clusive.  However,  for  plasmas  propagated  into  hydrogen  the  disagreement  between  the 
measured  conductivities  and  the  calculated  conductivities  was  very  striking.  In  the 
case  of  plasmas  extruded  into  a  3/8"  diamter  tube  the  results  were  o*  =0.13 
mho/cm.  o  Q  3  =  0.20  mho/cm,  Og^3  -  9.25  x  lO"^  mho/cm,  and  o^  ^3  =  0.11  x  lO'^  mho/cm 
giving  o  ^^3/0  Q  Q3  =  0.6  and  ^  -  82. 

(10)  Pressure  measurements  were  made  with  the  aquarium  technique  by  compressing 
plasmas  into  plastic  tubes  and  allowing  them  to  impact  water.  The  velocities  of  the 
transmitted  waves  in  water,  by  means  of  a  suitable  calibration  curve,  yielded  the 
pressure.  Measured  pressures  of  plasmas  were  found  in  this  manner  to  exceed  by  about 
a  factor  of  5  the  pressure  associated  with  a  shock  wave  considered  to  be  traveling 
at  the  velocity  of  the  plasma.  By  themselves  comparison  of  observed  and  computed 
pressure  probably  should  not  be  considered  sufficient  to  rule  out  the  possibility  of 
a  shock  mechanism.  However,  evidence  was  obtained  which  indicates  that  the  actual 
shock  wave  velocity  for  a  true  shock  wave  was  less  than  half  the  velocity  of  ionization 
waves  Hence,  the  actual  or  true  value  of  the  shock  wave  pressure  would  be  only  about 
1/4  as  great  as  the  value  computed  by  assuming  the  shock  wave  velocity  to  be  the  same 
as  the  ionization  wave.  On  this  basis  the  measured  ionization  wave  pressures  were 
10  to  20  times  greater  than  the  pressure  of  true  shock  waves  transmitted  from  the 
explosives  into  the  gaseous  medium.  The  velocity  of  the  observed  luminous  wave 
transmitted  from  Composition  B  directly  into  air  at  atmospheric  pressure  was  measured 
to  be  8.3  km/sec.  However,  when  an  aluminum  plate,  thin  enough  so  that  there  was  no 
appreciable  attenuation  as  determined  by  actual  pressure  measurements  by  the  aquarium 
method,  was  placed  on  the  end  of  the  charge  the  shock  wave  was  observed  to  travel 
at  only  4.4  km/sec.  Of  significance  is  the  fact  that  the  impedance  mismatch  equation 
was  found  to  predict  quite  accurately  the  shock  pressure  in  gaseous  media  of  such 
aluminum  terminated  charges  (even  though  he  equation  was  heretofore  not  beep  qon- 
sidered  applicable  for  such  large  differences  in  impedance).  Since  the  shock  impedance 
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of  aluminum  is  nearly  equal  to  that  of  Composition  B,  the  Impedance  mismatch  equation 
would  predict  very  nearly  the  same  shock  pressure  in  the  gas  for  bare  charges  as  for 
the  thin  plate  aluminum  terminated  charges. 

The  shapes  of  the  conduction  traces  measured  for  pure  shock  waves  were  different 
than  those  for  ionization  waves  from  bare  charges,  the  latter  exhibiting  a  more 
gradual  rise  at  the  front. 

(11)  Even  though  the  ionization  waves  from  hydrogen  and  helium  were  non- 
luminous  they  always  generated  high  luminosity  upon  impacting  the  end  of  the  glass 
or  plastic  chamber.  In  the  case  of  plasmas  from  Dithekite-13  (since  it  is  oxygen 
balanced  its  detonation  products  were  transparent)  close  observation  of  the  ionization 
waves  in  hydrogen  and  helium  revealed  a  diffuse,  non-luminous  zone  at  the  front  of 
the  detonation  products  that  corresponded  to  the  luminous  zone  in  air,  for  example. 
This  fuzzy  zone  was  believed  not  to  be  due  to  a  shock  wave  alone  since  it  was  too 
wide  to  be  attributed  to  possible  camera  "smear"  effects  and  it  did  not  exhibit  the 
change  in  index  of  refraction  associated  with  a  shock  wave.  The  fact  that  the 
fuzzy  zone  was  conducting  was  easily  determined  by  propagating  it  between  two  plates 
charged  to  a  potential  difference  of  20  KV  which  resulted  in  an  arc  discharge  through 
the  fuzzy  zone. 

It  was  concluded  finally  that  the  plasma  was  a  precursor  ionization  wave 
generated  primarily  from  chemionization  in  the  detonation  reaction  zone  which  is 
ejected  into  the  atmosphere  and  trapped  in  a  lower  density  state  when  the  electron- 
ion  recombination  times  are  longer.  However,  in  the  case  of  argon  or  air,  for 
example,  the  thermal  ionization  may  be  considerable.  The  ionized  material  is  then 
propagated  ahead  of  the  normal  shock  where  it  reacts  with  the  gaseous  medium  to 
produce  high  luminosity  in  many  cases.  It  is  not  unexpected  that  electrons  and 
perhaps  low  mass  ions  would  be  radiated  ahead  of  a  shock  owing  to  its  high  thermal 
gradient . 

The  luminosity  is  caused  by  the  interaction  of  electrons  in  the  plasma  primarily 
with  molecules  of  the  gaseous  medium  (largely  via  negative  ion  formation).  If 
negative  ions  are  not  obtained,  i.e.,  if  there  is  no  electron  affinity,  as  in  helium, 
no  medium  (vacuum)  or  if  there  is  rapid  ionization  decay  due  to  free  radicals,  there 
may  be  no  visible  or  only  a  transient  luminosity  associated  with  the  plasma  even 
though  it  actually  may  be  rich  in  free  electrons. 
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H.  The  Effect  of  Pressure  on  the  Degree  of  Ionization  in  Gaseous  Detonations 

The  results  of  this  study  were  presented  in  item  #11  of  the  list  of  technical 

notes.  In  recent  years  ionization  in  flames,  shock  waves,  and  in  detonation  waves 

has  been  the  subject  of  a  number  of  investigations.  The  Saha  equation 

log  k  »  log  [x/(l-x^)p]  -  -5024  E/T+  5/2  log  T+  log  g.g„/g„  -  6.47  (26) 

(x  is  the  fraction  of  the  species  ionized,  p  the  pressure,  E  the  ionization  potential, 

T  the  absolute  temperature,  and  g. ,  g  and  g  the  statistical  weights  of  the  ion, 

X  6  & 

electron,  and  atom,  respectively)  has  been  used  to  predict  equilibrium  constants 
for  reactions  of  the  type 

A  <=*  A*’  +  e‘ 

thus  permitting  calculations  of  electron  densities  in  thermal  ionization.  As  a 
general  rule  such  calculations  have  yielded  too  low  electron  densities  in  flames. 

For  detonations  in  condensed  explosives,  as  already  pointed  out,  the  calculated 
values  have  been  estimated  to  be  too  low  by  many  orders  of  magnitude.  On  the  other 
hand,  for  gaseous  detonations  in  hydrogen -oxygen  mixtures  at  relatively  low  pressure 
agreement  was  found  between  predictions  based  upon  Saha's  equation  and  experimental 
results,  although  at  50  atm  deviations  were  beginning  to  become  evident. 

A  large  and  Important  gap  in  experimental  electron  density  results  existed 
between  the  available  measurements  in  low  pressure  gaseous  detonations  (50  atm 
or  less)  and  results  for  condensed  explosives  where  detonation  pressures  are  up  to 
300  kb.  The  main  purpose  of  this  study  was  to  bridge  the  gap  as  completely  as 
possible  by  extending  electron  density  measurements  in  gaseous  detonations  to  as 
high  pressures  as  possible  and  to  determine  the  extent  of  departure  from  the 
thermodynamic  equilibrium  values  predicted  by  Saha's  equation. 

One  explanation  for  excess  ionization  as  pressure  is  Increased  was  based  upon 
the  quasi-lattice  model.  The  cohesive  energy  resulting  from  such  a  structure  should 
have  the  effect  of  lowering  the  effective  ionization  potential  of  the  species  in 
the  detonation  wave.  Thus,  as  pressure  or  density  increased  the  cohesive  energy 
would  Increase  and  the  ionization  potential  in  turn  would  decrease,  resulting  in  a 
greater  degree  of  ionization  than  predicted  by  Saha's  classical  equation.  Another 
explanation  attributed  the  excess  ionization  to  formation  of  hydrated  ions, 
principally  which  has  been  identified  by  mass  spectrometers  as  an  important 

ion  in  hydrogen-oxygen  flames. 
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Electron  density  measurements  were  determined  by  the  conduction  probe  method 

for  the  gaseous  mixtures  2H2  “  ~  0-002  N2,  2H2  “  ”  ^2  ”  ^2' 

pressure  range  covered  was  9.55  to  347  atm,  corresponding  to  initial  pressures  of 
0.54  to  17.35  atm  for  the  2H2  "  O2  -  0.002  N2  system,  17.3  atm  to  313  atm, 
corresponding  to  initial  pressures  of  0.99  atm  to  17.3  atm  for  the  2H2  "  ®2  "  ^2 
system,  and  for  2C0  -  0^  the  initial  pressure  was  varied  from  0.5  atm  to  10 
atm  which  corresponded  to  a  detonation  pressure  range  of  8.7  atm  to  191  atm. 

Thermohydrodynamic  calculations  were  carried  out  for  all  the  systems  correspond¬ 
ing  to  the  initial  pressures  of  interest.  Then  the  degree  of  ionization  of  each  of 
the  species  was  calculated  using  the  equilibrium  constants  determined  from  Saha's 
equation.  Due  to  the  low  ionization  potential  of  NO  it  was  essential  that  the 
concentration  of  N2  be  known  accurately  in  each  mixture.  In  the  mixture  2H2  "  ^2  " 
0.002  N2  the  0.002  N2  occurred  as  an  impurity  in  the  H2,  the  actual  concentration 
being  determined  by  analysis  with  a  mass  spectrometer.  For  the  hydrogen,  oxygen, 
nitrogen  mixtures  the  reactions  leading  to  ionization  which  were  first  considered 
were 

®2  °2^ 

H2O  H20'^  +  e’  (27) 

NO  NO"^  +  e” 

Ionization  contributed  by  the  other  products  of  detonation  H,  0,  H2  and  N2  was 
negligible  due  to  contributing  factors  of  high  ionization  potential  and  low  con¬ 
centration  of  the  species  in  the  detonation  products.  The  ionization  contributing 
reactions  considered  in  the  2C0  -  O2  system  were 

CO2  ^  C02^  +  e" 

CO  ^  CO"^  +  e’ 

(28) 

°2  °2 
0^0^+  € 

Comparisons  were  then  carried  out  between  the  measured  conductivities  and 
theoretical  conductivities  based  upon  the  calculated  electron  densities  according 
to  the  reactions  (27)  or  (28),  and  calculated  electron  mobilities.  In  the  case 
of  the  oxygen-hydrogen-nitrogen  mixtures,  the  theoretical  and  experimental  con¬ 
ductivities  were  found  to  agree  at  the  lowest  pressures  for  which  data  were  obtained. 
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For  the  2CO  -  0^  system,  however,  the  experimental  conductivities  exceeded  the 

calculated  conductivities  by  roughly  50  times.  Thus,  apparently  some  important 

ion  contributing  species  must  have  been  neglected  in  this  system.  For  both  the 

hydrogen-oxygen-nitrogen  mixtures,  as  pressure  was  increased,  the  electrical 

conductivity  was  found  ‘o  rise  sharply  above  the  theoretical  curve  calculated  on 
+  -I- 

the  basis  of  H2O  ,  O2  and  NO  formation,  indicating  therefore  that  an  additional 
mechanism  contributing  ionization  became  Important.  Studies  with  mass  spectro¬ 
meters  had  revealed  H_0^  to  be  an  important  ion  species  in  hydrogen-oxygen  flames, 

+ 

and  since  H^O  formation  would  be  expected  to  become  more  pronounced  at  high  pressures 
theoretical  analyses  were  carried  out  to  determine  the  influence  upon  electrical 
conductivity  when  was  considered. 

An  equilibrium  constant  for  the  reaction 

H20'^  +  H  (29) 

was  calculated  from  partition  functions  using  an  enthalpy  of  reaction, at 
298°K  equal  to  -90  kcal/mole  which  was  considered  to  be  the  best  value  available. 

This  equilibrium  constant  turned  out  to  be 


0.197  X  10‘^°  T"^/^  ^44,840/T  _  /(H  (H)  (molecules/cm^) 


(30) 


^37  . .  ‘  '•**3''  "'>•‘2' 

Reaction  (29)  was  then  considered  simultaneously  with  reactions  (27)  and  new  electron 

densities  and  corresponding  conductivities  were  calculated.  Although  the  contribution 

to  the  total  ionization  by  the  formation  of  H^o"*”  Increased  with  pressure,  this 

contribution  was  much  less  than  was  required  for  agreement  with  experimental  results. 

For  example,  at  about  250  atmospheres  using  the  2H-  -  0-  -  0.002  N  system 
+  2  2  2 

inclusion  of  H^O  in  the  theoretical  calculations  Increased  the  degree  of  ionization 
only  a  negligible  amount  whereas  approximately  a  tenfold  increase  was  required. 

In  order  to  determine  if  the  inclusion  of  H^o"*”  formation  in  the  equilibria 
could  possibly  account  for  the  sharp  increase  in  the  degree  of  ionization  observed 
as  pressure  was  Increased  above  20  atmospheres,  a  value  ofAH°  for  the  reaction 
given  by  equation  (29)  was  chosen  such  that  the  calculated  and  experimental  values 
of  the  conductivity  agreed  at  520  atm.  This  value  of  AH°  turned  out  to  be 
-172.75  kcal/mole  and  the  corresponding  equilibrium  constant  was 


=  0-197  X  10‘^°  .j,-3/2  ^  86,300/T 


(31) 


In  this  instance  the  calculations  predicted  much  too  high  a  degree  of  ionization 
for  the  low  pressure  region.  Aside  from  the  fact  that  AH®  =  -172.72  kcal/mole 
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certainly  was  an  unacceptable  value,  the  mechanism  of  H^O  formation  also  proved 
unable  to  produce  the  observed  shape  of  the  electrical  conductivity  vs  pressure 
curve . 

The  quasi-lattice  model  was  then  considered  as  a  possible  explanation  of  the 
excess  ionization.  In  this  model,  provided  the  electron  density  is  sufficiently 
large,  the  ions  and  electrons  are  considered  to  assume  a  quasl-metalllc-llke  lattice 
structure  in  which  the  electrons  move  through  the  positive  ion  lattice  in  correlated 
patterns.  Such  a  situation  would  lead  to  a  cohesive  energy,  E^,  given  approximately 
by  the  equation  for  metals 

4/5  o 

-E  =  BZ(l-0.8/dz  '•^)/d  (e.v.  for  d  in  A)  (32) 

P 

where  d  is  the  lattice  distance,  Z  is  the  number  of  electrons  per  ion,  and  B  is 
an  emperlcal  constant.  This  equation  would,  of  course,  apply  only  for  situations 
where  -E^  was  sufficiently  greater  than  kT  that  the  quasi-lattice  structure  was 
s  table . 

Due  to  the  existence  of  cohesion,  Saha's  equation  for  thermionlzatlon  was 
modified  by  subtracting  from  the  E  term  the  cohesive  energy  E^ ,  and  calculations 
of  the  degree  of  ionization  were  carried  out  for  the  2H2  "  ”  ^2  with 

the  new,  resulting  equilibrium  constants  using  various  values  of  the  emperlcal 
constant,  B.  Using  the  theroetical  value  of  B  equal  to  30,  which  had  been 
determined  for  metals,  yielded  a  conductivity  vs  pressure  curve  that  possessed 
the  same  shape  as  the  experimental  curve,  but  the  electron  densities  predicted 
were  too  high.  However,  use  of  a  value  of  B  about  1/3  of  this  value  gave  good 
agreement  throughout  the  whole  range  of  pressures. 

The  sharp  increase  in  conductivity  with  pressure  observed  in  the  hydrogen, 
oxygen,  nitrogen  systems  was  not  found  with  2C0  -  O2  contrary  to  expectations 
based  upon  the  quasi- lattice  model.  However,  in  this  case  theoretical  and  experimental 
results  proved  to  deviate  widely  even  at  low  pressure.  Thus,  one  did  not  have  a 
reliable  theoretical  reference  for  low  pressure  from  which  to  judge  how  theory  and 
experiment  deviated  at  high  pressure. 

I ,  Ionization  and  Electrical  Conductivity  and  its  Relationship  to  the 
Deflagration  to  Detonation  Transition  in  Solid  Explosives 

At  the  preparing  of  this  final  report  the  investigations  concerned  with  the 
above  program  are  just  being  terminated  and  a  technical  note  covering  the  study 
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wlll  be  prepared  in  the  near  future.  Hence,  only  a  description  of  the  problem, 
its  background,  and  a  few  of  the  results  will  be  elucidated  here. 

The  classical  theory  of  detonation  is  based  upon  transport  phenomena  being 
negligible  throughout  the  reaction  zone.  According  to  this  theory  a  detonation 
wave  is  simply  a  shock  wave  supported  by  chemical  reaction,  the  chemical  reaction 
being  initiated  at  the  wave  front  by  the  Hugionot  heat  of  compression  of  the  shock 
wave.  This  theory  predicts  a  pressure  "spike"  at  the  wave  front  whose  magnitude 
is  approximately  twice  the  Chapman  Jouguet  pressure,  which  in  the  case  of  ideal 
explosives,  is  the  pressure  at  the  end  of  the  reaction  zone.  The  pressure  at  the 
wave  front  decays  through  the  reaction  zone  at  a  rate  proportional  to  the  chemical 
reaction  rate  along  the  so-called  Todes  or  Raleigh  line. 

Measurements  of  pressures  at  detonation  wave  fronts  performed  by  this  group 
under  the  auspices  of  another  contract^ negated  the  presence  of  a  spike,  however. 
The  pressure  at  the  detonation  wave  front  was  found  not  to  exceed  the  Chapman 
Jouguet  pressure.  These  determinations  incidentally  were  performed  for  explosives 
possessing  a  wide  variety  of  reaction  zone  lengths,  the  longer  reaction  zone  lengths 
being  sufficient  that  according  to  any  existing  theory  a  "spike"  should  have  easily 
been  discernable.  Studies  under  still  another  contract'  indicated,  in  fact,  that 
the  pressure  exhibited  a  finite  rise  time  at  the  front.  According  to  the  classical 
model,  the  temperature  at  the  wave  front  would  only  be  the  order  of  several  hundred 
degrees  and  thus  the  chemical  reaction  should  be  slow  at  first.  This  fact  it  would 
seem  should  lead  to  excessively  long  reaction  times  and  a  chemical  induction  period 
at  the  wave  front. 

The  presence  of  the  high  electron  densities  in  the  detonation  reaction  zones 
as  described  earlier  in  the  report,  however,  provided  a  mechanism  whereby  conduction 
of  heat  in  the  reaction  zone  could  be  considerable  owing  to  the  high  mobility  of 
electrons  under  the  influence  of  a  sharp  thermal  gradient  at  the  wave  front.  Thus, 
conceivably  appreciable  amounts  of  heat  could  be  conducted  from  within  the  reaction 
zone  to  the  wave  front,  the  effects  of  which  would  be  to  diminish  magnitude  of  the 
"spike"  erasing  it  completely  if  the  heat  conduction  were  sufficient,  and  to  increase 
the  temperature  at  the  detonation  wave  front  thereby  increasing  the  chemical  reaction 
rate  at  the  front. 

High  electron  densities  and  apparently  high  heat  conduction  therefore  were 
existing  conditions  in  detonation.  The  next  question  of  interest  was,  naturally. 
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are  they  necessary  conditions  for  detonation?  This  question  prompted  the  initiation 

of  the  program  with  which  this  section  is  concerned.  As  a  part  of  another 

investigation  concerned  with  the  transition  to  detonation  for  the  shock  initiation 

of  liquid  explosives  in  the  so-called  "card-gap"  test,  photographs  showed  a  super 

velocity  phenomenon  of  low  luminosiry  starting  at  the  barrier-receptor  Interface 

and  terminating  at  the  initial  sh..>ck  front  with  the  transition  to  high  order 

(81 

detonation  occurring  at  the  intersection  of  these  two  phenomena.^  ■'  In  line  with 
the  ideas  of  the  importance  of  ionization  in  the  detonation  reaction  zone  this 
super  velocity  phemomer.a  was  interpreted  to  be  a  wave  of  ionization  which  cascaded 
to  the  shock  front  from  the  barrier  receptor  interface  as  a  result  of  chemical 
reaction  initiated  by  the  initial  shock  wave  and  was  termed  the  "flash  across"  or 
"heat  pulse". 

The  studies  with  which  this  section  pertains  involved  the  measurements  of 
conductivity  for  solid  explosives  initiated  by  shock  in  the  "card  gap"  test  and 
Interpreting  the  relationship  between  conductivity  and  the  deflagration  to  deto¬ 
nation  transition.  The  measurements  were  made  with  the  conduction  probe  method 
primarily  in  the  region  between  the  barrier-receptor  interface  and  the  distance, 

S2»  beyond  it  where  high  order  detonation  first  occurred,  both  single  sets  and 
double  sets  of  perpendicular  probes  as  well  as  parallel  probes  being  employed. 

The  perpendicular  probes  were  so-called  because  the  conductors  were  perpendicular 
to  the  charge  axis  with  the  probe  gap  on  the  charge  axis  thus  providing  a  measure¬ 
ment  of  conductivity  vs  time  at  one  point  within  the  charge  for  each  set  of 
probes.  The  parallel  probe  system,  on  the  other  hand,  consisted  of  two  conductors 
mounted  symmetrically  about  and  parallel  to  the  charge  axis,  thus  giving  an 
integrated  value  of  the  conductivity  between  them.  The  main  use  of  the  double 
probe  system  was  to  verify  in  which  direction  (forward  or  backward)  each  of  the 
disturbances  propagated. 

The  investigations  were  commenced  under  the  contract  that  this  final  report 
covers,  namely  Contract  AF-18(603)-100.  However,  they  continued  into  the  present 
contJact,  AF-49(638) 1061 ,  and  at  present  have  reached  the  logical  state  where  a 
technical  note  should  be  prepared  soon.  One  of  the  most  important  results  that 
was  determined  in  this  program  was  confirmation  of  the  "flash  across"  mechanism. 

In  the  transition  state  before  high  order  detonation  has  been  established  the 
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ionization  wave  was  found  to  exhibit  a  long,  gradual  rise.  The  ionization  started 
apparently  at  the  pressure  front  but  during  the  initial  stages  the  rise  time  of  the 
ionization  wave  was  much  longer  than  the  rise  time  of  the  initial  pressure  wave. 
However,  as  the  disturbance  propagated  through  the  donor  charge  the  ionization  wave 
was  found  to  sharpen,  and  after  the  ionization  wave  had  sharpened  to  a  steep  fronted 
wave  which  coincided  apparently  with  the  Initial  pressure  wave  the  transition  to 
high  order  detonation  occurred.  In  other  words,  the  transition  took  place  when  the 
ionization  front  merged  with  the  shock  front  permitting  conduction  of  heat  to  the 
shock  front  One  would  expect  that  a  critical  amount  of  ionization  would  be  required 
before  the  associated  luminosity  would  be  sufficient  for  photography.  Thus,  the 
super  velocity  "flash  across"  phenomenon  observed  in  liquids  probably  corresponded 
to  the  propagation  of  a  critical  level  of  ionization.  The  velocity  associated  with 
it  would  be  considerably  greater  than  the  velocity  of  the  initial  pressure  wave  as 
a  result  of  the  "sharpening  up"  of  the  ionization  wave. 


•  Transition  to  Detonation  in  Liquid  Explosives 
As  mentioned  in  the  previous  section  a  super  velocity,  slightly  luminous 
phemomenon  has  been  observed  to  start  at  the  barrier-receptor  interface  and  propagate 
to  the  initial  shock  front  at  which  point  normal  high  order  detonation  is  established. 


The  interpretation  advanced  at  this  laboratory  was  that  the  phemomenon  was  not  of 
a  "hydrodynamic"  nature  but  corresponded  to  a  critical  level  of  ionization  which  by 
virtue  of  the  chemical  reaction  started  by  the  initial  shock  and  by  heat  conduction, 
cascaded  to  the  initial  shock  front,  high  order  detonation  occurring  when  the 
ionization  wave  with  its  associated  high  conductivity  merged  with  the  initial  shock. 
Another  explanation  has  been  advanced  for  this  phenomenon, namely  that  the 
super  velocity  phenomenon  is  a  high  order  detonation  which  starts  at  the  barrier- 
receptor  interface  propagating  at  super  detonation  velocity  to  the  shock  front  by 


virtue  of  the  high  density  of  the  compressed  medium  between  the  initial  shock  front 
and  the  barrier.  This  explanation  required  that  the  mechanism  of  shock  initiation 
is  different  in  liquid  than  in  solids  for  which  the  transition  to  detonation  occurs 
some  distance  within  the  receptor  A  crucial  experiment  which  would  determine 
whether  or  not  the  latter  explanation  is  correct  is  the  measurement  of  the  pressure 
associated  with  the  super  velocity  phenomenon.  If  it  were  indeed  a  detonation  the 
pressure  would  be  exceedingly  high,  having  a  value  several  times  the  normal 
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detonation  pressure  of  the  liquid,  whereas  if  it  were  an  ionization  phenomenon, 
the  pressure  would  be  lower  than  the  pressure  of  the  initial  shock. 

The  task  of  measuring  the  pressure  associated  with  the  t.ypervelocity  phenomenon 
in  liquids  was  undertaken  during  the  latter  part  of  the  contract  period,  and  the 
investigations  are  being  continued  under  the  auspices  of  the  new  contract.  The 
only  method  of  estimating  the  pressure  which  would  appear  to  give  valid  results  is 
the  aquarium  technique  in  which  the  shock  is  transmitted  into  an  inert  medium  of 
known  equation  of  state.  The  use  of  such  a  method  requires  especially  careful 
techniques  in  this  case.  The  transition  to  detonation  or  the  overtaking  of  the 
initial  shock  by  the  detonation  wave  that  started  at  the  plate,  whichever  the  case 
may  be,  should  occur  slightly  inside  pressure  measuring  medium  (imagining  that  the 
pressure  measuring  medium  were  the  same  as  the  explosive).  Hence,  it  is  important 
that  the  shock  Impedance  of  the  explosive  match  the  impedance  of  the  transmitting 
medium  in  order  that  reflections  do  not  occur  at  the  interface  which  might  influence 
the  transition.  Also,  is  necessary  that  distance  from  the  barrier  at  which  the 
coelescence  occurs  be  known  accurately  so  that  the  proper  length  of  receptor  charge 
be  known.  If,  then,  the  hypervelocity  phenomenon  is  a  detonation,  a  strong  shock 
velocity  "overshoot"  must  be  measured  in  the  transmitting  medium. 

At  the  end  of  the  contract  period  positive  results  had  not  been  obtained. 
However,  the  Sj^  vs  $2  relationships,  where  Sj^  is  the  barrier  thickness  and  S2  is 
the  distance  beyond  the  barrier  at  which  the  hypervelocity  phenomenon  and  the 
initial  shock  front  merged,  had  been  determined  for  nitromethane  using  a  particular 
charge  geometry,  and  a  transmitting  medium  whose  equation  of  state  was  a  good  match 
for  MM  has  been  found.  During  the  new  contract  period  efforts  will  be  devoted  to 
this  problem  through  studies  with  a  number  of  liquids.  Determination  of  the  correct 
interpretation  of  the  "flash  across"  phenomenon  would  be  an  exceedingly  important 
contribution  to  the  theory  of  detonation. 
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